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ABSTRACT 

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi- 
wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope 
mission. While investigating the previously reported 7 -ray/optical flare accompanied by a change in 
optical polarization, we found that the optical emission appears delayed with respect to the 7 -ray 
emission by about 10 days. X-ray observations reveal a pair of ‘isolated’ flares separated by ~ 90 
days, with only weak 7 -ray/optical counterparts. The spectral structure measured by Spitzer reveals 
a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared baud 
during the 7 -ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected 
spectral energy distributions are fitted with leptonic models including Comptonization of external 
radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the po- 
larization swing involving propagation of the emitting region along a curved trajectory, we can explain 
the evolution of the broad-band spectra during the 7 -ray flaring event by a shift of its location from 
~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the 7 - ray flare is generated instead 
at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron 
self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking syn- 
chrotron component, and suggest that the corresponding inverse- Compton component explains the 
steady X-ray emission. 

Keywords: radiation mechanisms: non-thermal — galaxies: active — galaxies: jets — quasars: indi- 
vidual (3C 279) — gamma rays: galaxies — X-rays: galaxies 
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1. INTRODUCTION 

Blazers are Active Galactic Nuclei (AGN) character- 
ized by highly luminous and rapidly variable continuum 
emission at all observed bands. The most commonly ac- 
cepted scenario has their broad-band emission Doppler- 
boosted by a relativistic jet pointing close to our line of 
sight (e.g., Ulrich et al. 1997). While the jet emission 
usually dominates the observed broad-band spectrum, 
the optical/ultra-violet (UV) and infrared (IR) spectra 
often also reveal signatures of the central engine: broad 
emission lines, and in some cases, quasi-thermal opti- 
cal/UV emission and IR dust emission, indicating the 
presence of an accreting supermassive black hole. Most 
viable current models for the origin of such jets involve 
conversion of the gravitational energy of matter flowing 
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onto the black hole to the kinetic energy of the relativistic 
outflow or tapping the rotation energy of a spinning black 
hole. However, the conversion process itself is not well 
understood, and many additional questions regarding the 
dissipation region of the jet’s energy into radiation and, 
in particular, its location remain unanswered. 

Major advances in understanding of blazars came as 
a result of the discovery by the EGRET instrument on 
board Compton Gamma- Ray Observatory (CGRO) that 
they are strong 7-ray emitters, with 7-rays dominat- 
ing radiative output (Hartman et al. 1999). With this, 
multi-band observations including the 7-ray band hold 
the promise of answering many outstanding questions re- 
garding the structure of the relativistic jets of blazars. 

3C 279 ( z = 0.536; Lynds et al. 1965) is in fact 
one of the first 7-ray blazars discovered by EGRET in 
1991 (Hartman et al. 1992). The 7 -ray signal had been 
significantly detected in each observation by EGRET 
since its discovery (see, e.g., Hartman et al. 2001a), with 
the flux having ranged over roughly 2 orders of magni- 
tude, from ~ 10 -7 up to ~ 10 -5 photons cm -2 s -1 above 
100 MeV (Maraschi et al. 1994; Wehrle et al. 1998), and 
a factor of 2 variation on timescales as short as 8 hrs. 
The photon index in the EGRET 7-ray band ranged 
from 1.8 to 2.3 (Nandikotkur et al. 2007). On a few oc- 
casions, 3C 279 was also detected at lower energies by 
G’GRO’.sOSSE (50 keV-1 MeV) (McNaron-Brown et al. 
1995) and COMPTEL (0.75-30 MeV) (Hermsen et al. 
1993; Collmar et al. 2001) instruments, indicating that 
the 7 -ray emission forms a broad peak in the vF v rep- 
resentation. In 2008 July, the AGILE satellite ob- 
served a 7-ray flare associated with the source with 
11.1 p significance (Giuliani et al. 2009), with an av- 
erage flux above 100 MeV of (21.0 ± 3.8) x 10 -7 
photons cm -2 s -1 and the photon index of 2.22±0.23 be- 
tween 100 and 1000 MeV. In the very-high-energy (VHE) 
7-ray regime above 100 GeV, the imaging atmospheric 
Cherenkov telescope MAGIC detected flares twice in 
2006 February (Albert et al. 2009) and in 2007 Jan- 
uary (Aleksic et al. 2011), which made this source the 
most distant currently known VHE 7-ray emitter. 

Optical and UV observations of the source in relatively 
low states - when the jet emission was relatively faint - 
allowed a study of the accreting black hole and the associ- 
ated accretion disk. The luminosity of the accretion disk 
was estimated to be Ld ~ 2 x ICr 8 erg s -1 (Pian et al. 
1999). The mass of central supermassive black hole was 
estimated to be in the range of (3-8 )x10 8 Mq using the 
luminosity of optical broad line (Woo & Urry 2002) or 
the Hg line width (Gu et al. 2001). Those values are 
similar to the estimates based on the luminosity of the 
host galaxy (Nilsson et al. 2009). 

3C 279 contains a compact radio core, associated with 
time-variable jet-like structure. Radio observations at 
43 GHz by Very Long Baseline Array (VLBA) between 
1998 March and 2001 April revealed superluminal mo- 
tion of the jet with apparent speeds that range from 5c 
to 17c (Jorstad et al. 2004, 2005). Those observations 
also allowed an estimate of the Lorentz factor of the jet 
flow of Tj = 15.5 ±2.5 and of the viewing angle of the jet 
0 O = 2.1 ± 1.1 deg, corresponding to a Doppler beaming 
factor of 6 = 24.1 ± 6.5. A change of the trajectory of 
a jet component has also been reported in radio obser- 
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vations with VLB A (Homan et al. 2003). Those authors 
estimated the jet component to be moving with a Lorentz 
factor Tj > 15 at an initial viewing angle of < 1° . 

The broad-band spectral energy distribution (SED) 
of the source is characterized by a two-bump struc- 
ture, similar to many other 7 - ray blazars. In the 
context of widely accepted leptonic models, the lower- 
frequency bump, peaking at the far-IR and extend- 
ing to the extreme UV band, is commonly ascribed to 
synchrotron radiation from relativistic electrons in the 
jet. The second bump, spanning from the X-ray to 
the 7 -ray band with a peak in the MeV-GeV range, 
is believed to be generated via inverse- Compton scat- 
tering, presumably by the same population of parti- 
cles that radiate at lower energies via the synchrotron 
process. The seed photons for the Compton scattering 
can be synchrotron photons (synchrotron self-Compton: 
SSC, Maraschi ct al. 1992; Bloom & Marscher 1996), 
accretion disk photons (external Compton scattering 
of direct disk radiation: ECD, Dermer et al. 1992; 

Dermer & Schlickeiser 1993) and accretion disk photons 
re-scattered by the broad-line region clouds/intercloud 
meciium (ECC, Sikora et al. 1994; Blandford & Levinson 
1995), or infrared radiation from a torus located be- 
yond the broad- line region (ERC-IR, Sikora et al. 1994). 
Specific to 3C 279, multiwavelength snapshot observa- 
tions for severed epochs including 7 -rays were presented 
in Hartman et al. (2001a). Those authors explained the 
overall spectra using the leptonic model, where the X-ray 
photons are mainly produced by SSC, and both ECD and 
ECC contribute to the 7 -ray emission. Spectral variabil- 
ity was explained by variations of the bulk Lorentz factor 
of the jet, accompanied by changes in the spectral shape 
of the electron distribution. 

The optical variability of 3C 279 is extreme: in 1937, it 
showed optical ( B ) magnitude of 11.27 (Eachus & Liller 
1975), making it one of the most luminous active galax- 
ies ever recorded. The strong variability recorded in all 
bands provides an opportunity to establish the relation- 
ship between emission in those bands, and thus can be 
used to constrain theoretical models of physical regions of 
the jet responsible for such emission. Many such multi- 
wavelength campaigns have been conducted (see, e.g., 
Maraschi et al. 1994; Wehrle et al. 1998; Larionov et al. 
2008; Collmar et al. 2010) but they have not revealed 
a simple relationship between the variability in various 
bands: radiation in different spectral regimes does not 
always rise and fall simultaneously, although the periods 
of increased rapid activity in all bands seem to last for 
several months, and take place when the source is rela- 
tively bright. A recent paper by Chatterjee et al. (2008) 
presents the results of the monitoring of 3C 279 for 11 
years in radio, optical and X-rays, and discusses the de- 
tails of the jet structure based on multi-band correlation 
studies. However, due to the lack of deployed instru- 
ments, long term monitoring observations could not in- 
clude the 7 -ray regime, where the source often shows 
stronger variability than in other bands. 

The launch of the Fermi Gamma-ray Space Telescope 
on 2008 June 11 has rejuvenated multi-band studies of 
blazars The Large Area Telescope (LAT; Atwood et al. 
2009) instrument on Fermi can monitor all 7 -ray sources 
on the sky with its a wide filed of view and a much larger 
effective area compared to earlier 7 -ray missions. Taking 


advantage of this new instrument for 7 -ray observations, 
we have organized intensive multiwavelength campaigns 
for 3C 279 from radio to the high-energy 7 -ray energy 
ranges. Many ground-based telescopes (cm, mm, near- 
IR and optical) and various satellites (IR, UV, X-ray, 
hard X-ray and high-energy 7 -ray) participated in this 
campaign. We reported the first results of the campaign 
in Abdo et al. (2010a, hereafter Paper I), where we dis- 
covered the dramatic change of the optical polarization 
coincident with the 7 -ray flare. Here, we provide details 
of the multi-band observations and the interpretation of 
those data for the 2-year interval between 2008 August 
and 2010 August. In Section 2, we present and briefly 
discuss the features of the LAT 7 -ray data; in Section 
3, we present the data in lower energy bands. Section 4 
highlights the features of time series measured in various 
bands including their cross correlations, and the general 
properties of the broad-band spectral energy distribu- 
tion. In Section 5, we provide viable emission models for 
the source in the context of leptonic scenarios. 

2, Fermi-LAT DATA AND RESULTS 

Fermi-LAT is a pair-production telescope with large 
effective area (8,000 cm 2 on axis at lGeV for the event 
class considered here), and large field of view (2.4 sr at 
lGeV), sensitive to 7 rays in the energy range from 
20 MeV to > 300 GeV. Information regarding on-orbit 
calibration procedures is given in Abdo et al. (2009a). 
Fermi-LAT normally operates in a scanning ‘sky-survey’ 
mode, which provides a full-sky coverage every two orbits 
(3 hours). For operational reasons, the standard rocking 
angle (defined as the angle between the zenith and the 
center of the LAT field of view) for survey mode was 
increased from 35° to 50° on 2009 September 3. 

2 . 1 . Observation and data reductions 

The data used here comprise 2-year observations ob- 
tained between 2008 August 4 and 2010 August 6 (MJD 
54682 - 55414). We used the standard LAT analysis soft- 
ware, ScienceTools v9r21. The events were selected using 
so-called “diffuse class” events. In addition, we excluded 
the events with zenith angles greater than 100 ° to avoid 
the contamination of the Earth-limb secondary 7 radi- 
ation. The events were extracted in the range between 
200 MeV and 300 GeV within a 15° acceptance cone of 
the Region of Interest (ROI) centered on the location of 
3C 279 (RA = 195.047°, DEC=-5.789°, J2000). Below 
200 MeV, the effective collection area of LAT for the dif- 
fuse class events drops very quickly and thus larger sys- 
tematic errors are expected. The 7 -ray flux and spectrum 
were calculated using the instrument response function 
(IRF) of “P 6 _V 11 -DIFFUSE” by an unbinned maximum 
likelihood fit of model parameters. We examined the sig- 
nificance of the 7 -ray signal from the sources by means 
of the test statistic ( TS ) based on the likelihood ratio 
test 59 . The background models included a component 
for the Galactic diffuse emission along the plane of the 
Milky Way, which was modeled by the map cube file 
“ glLiem-vO$LP6- V1 1 .DIFFU SE.fitd ' . An isotropic com- 
ponent ( isotropic-iem-v02-P6-Vll-DIFFUSE.txt ) was 

59 TS = 25 with 2 degrees of freedom corresponds to an es- 
timated ~ 4.6 ct pre-trials statistical significance assuming that 
the null-hypothesis TS distribution follows a x 2 distribution (see 
Mattox et al. 1996). 
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also included to represent the extragalactic diffuse emis- 
sion and residual instrumental background. Besides 
those components, the model in our analysis also in- 
cluded the emission from all nearby point sources in- 
side the ROI from the first Fermi-LAT catalog (1FGL: 
Abdo et al. 2010b). The spectra of those sources were 
modeled by power law functions except for a pulsar 
1FGL J 1231. 1-1410 (=PSR J1231-1411), for which we 
included an additional exponential cut-off in its spectral 
modeling. During the spectral fitting, the normalization 
factors of the Galactic diffuse and isotropic components 
and the nearby sources were left as free parameters, and 
the photon indices of the nearby sources were fixed to the 
values from the 1FGL catalog except for 3C 273, whose 
photon index was allowed to vary freely. In the light 
curve analysis, we considered only two bright sources in 
the background model as nearby point sources, namely 
3C 273 and 1FGL J1231. 1-1410, because other nearby 
sources had a negligible contribution to 7-ray signal, es- 
pecially in such relatively short time scales (shorter than 
a week) for the light curves considered here. The fluxes 
used for the light curve were calculated by a simple power 
law model fit using data in the given energy ranges. 


2.2. Temporal behavior 

The 7-ray light curve measured by Fermi - LAT can be 
seen in Figure 1. The Figure shows the flux history above 
200 MeV averaged over (a) 1-day intervals, (b) 3-day in- 
tervals, and (c) 1-week intervals. It also includes 1-week 
light curves of (d) the flux between 200 MeV and 1 GeV, 
(e) the flux above 1 GeV, and (f) the photon index in 
the range above 200 MeV. 

The 7-ray flux clearly shows variability. The source 
showed high-flux states between MJD 54700 and 54900, 
in which two prominent flares can be seen: one of 
the flares at ~ MJD 54800 and the the other at ~ 
MJD 54880. During the second flare, a change in the 
optical polarization associated with a 7-ray flare was 
discovered (Paper I). We detected some flux variabil- 
ity between MJD 55000 and 55120, but after that, the 
source remained in a relatively low-activity state until 
the end of the period considered in this paper. Dur- 
ing this 2-year period, the highest integral flux above 
200 MeV occured oh MJD 54880 in the 1-day interval 
light curve with flux of Ffi>200MeV= (11.8 ± 1.5) x 10 -7 
photons cm -2 s _1 and TS = 306. By extrapolating 
the spectrum down to 100 MeV, an integral flux above 
100 MeV on that day yields Fe>iooMoV= (31.0 ± 6.0) x 
10~ 7 photons cm -2 s -1 , v/hich is still a factor of 3-4 
times lower than the flux of the brightest flare (~ 1 x 10 -5 
photons cm -2 s _1 ) detected during the EGRET obser- 
vations of the source (Wehrle et al. 1998; Hartman et al. 
2001b). 

We quantified the flux variability using 1-week inter- 
val data for energies above 200 MeV [full band], between 
200 MeV and 1 GeV [soft band] , and above 1 GeV [hard 
band]. This is based on the “excess variance'’ method 
(Nandra et al. 1997; Edelson et al. 2002) after subtract- 
ing the contribution expected from measurement errors 
(uerr,*)- Using the mean square error < <7 e rr,i >1 the ex- 
cess variance Evar can be described as (Vaughan et al. 


2003) 




S 2 — < U e 




<F > 2 


( 1 ) 


where S is the variance of the flux, and < F > is the 
mean value of the flux. The definition of associated error 
can be found in (Vaughan et al. 2003). In the calcula- 
tion, we excluded bins of 8, 82, 87 and 90 because the 
fit in the hard band failed due to poor statistics of the 
data samples. Resulting F w values are 0.695 ± 0.015, 
0.648 ±0.017 and 0.839 ±0.030 for the full, soft and hard 
bands, respective^*. The resulting values indicate that 
the flux of the hard band showed significantly stronger 
variability than that of the soft band. For comparison, 
F, a r = 0.79 ± 0.02 for E > 300 MeV has been reported 
during the first 11 months of the Fermi scientific mis- 
sion (Abdo et al, 2010f), when the source has clearly 
been more active. 

A Power Density Spectrum (PDS) for the 3-day binned 
light curve was calculated using a Fourier transform 
and is shown in Figure 2. The power density was 
normalized to fractional variance per frequency unit 
( rms 2 I~ 2 Day -1 ) and the PDS points were averaged in 
logarithmic frequency bins. The white noise level was 
estimated from the rms of the flux errors and was sub- 
tracted from the PDS. A slope of 1.6 ± 0.2 was obtained 
from a linear fit to the binned PDS for frequencies up 
to 0.1 Day -1 . The main uncertainty in the estimated 
PDS slope is due to the stochastic nature of the variabil- 
ity which leads to variations in the determined slope be- 
tween different time limited observations. An additional 
effect which can cause a systematic bias in the observed 
PDS slope is the red noise leakage (e.g., Chatterjee et al. 
2008). In the present analysis this effect is not taken into 
account. 

Figure 3 shows plots of flux vs. photon index (T) 
based on the weekly light curve results above 200 MeV 
[full band], between 200 MeV and lGeV [soft band] 
and above 1 GeV [hard band] . The data which have 
TS > 10 were selected for the plots and are shown in 
gray points. An average photon index was calculated by 
fitting a constant value in each plot, corresponding to 
I* >200 MeV = 2.334 ±0.015, F200MeV-lGeV = 2.20 ± 0.03 
and r> i GeV = 2.48 ± 0.04 for the full, soft and hard 
bands, respectively. The average photon index in the soft 
band shows a significantly harder spectrum than that in 
the hard band. 

We also derived photon indices resulting from an 
analysis where the data were sorted in five bins using 
week-long fluxes for each energy band, and plotted the 
results as red points. Those photon indices of each flux 
bin are also shown in the insets in Figure 3. For the 
full band, although the change of the photon index is 
rather small (Ar ~ 0.2) compared to the flux variation 
(spanning about an order of magnitude), a weak “harder 
when brighter” effect can be seen. Such an effect 
was also measured in other LAT blazars (Abdo et al. 
2010d). The soft band also shows the weak “harder 
when brighter” effect with a slightly larger change of 
the photon index (Ar ~ 0.4). On the other hand, the 
photon index of the hard band changes only sightly 
(Ar ~ 0.1) and is statistically consistent with a constant 
value. 
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Figure 1. Gamma-ray light curves of 3C 279 during the first two years of the Fermi- L AT observations from 2008 August to 2010 August, 
plotted in; (a) 1-day intervals at energies above 200 MeV, (b) 3-day intervals at energies above 200 MeV, (c) 1-week intervals at energies 
above 200 MeV, (d) 1-week intervals at energies between 200 MeV and 1 GoV, (e) 1-week intervals at energies above 1 GeV. The panel-(f) 
shows the history of the photon index at energies above 200 MeV in 1 week intervals, while the panel-(g) shows arrival time distribution 
of > 20 GeV events associated with 3C 279. The vertical axis of the panel-(g) represents the estimated energy of events. The highest 
energy photon corresponds to 30.8 GeV at MJD 54891. The dotted lines and capital letters represent time intervals where 7 -ray spectra 
are extracted (see also Table 1). 
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Figure 2. Power Density Spectrum of 3C 279 for the 3-day binned 
7 -ray light curve. The white noise level has been subtracted. The 
solid line histogram describes the PDS averaged in logarithmic fre- 
quency bins while the dotted curve describes the raw PDS before 
binning. The dashed line represents a linear fit to the binned PDS. 


2.3. Highest energy photons 

During the 2-year observations, the highest en- 
ergy photon associated with 3C 279 was detected at 
MJD 54891.60745 with an estimated energy of 30.8 GeV. 
The event was converted in the front-thin layers (so- 
called “front event”) of the LAT detector and still re- 
mains even when we apply the cleanest event selec- 
tion (so-called “data clean event”), which was developed 
for studying extragalactic 7 -ray background (Abdo et al. 
201 Og). The reconstructed arrival direction of the event 
is 5' .7 (=0°.095) away from 3C 279, and is within the 68 % 
containment radius of the LAT PSF (0°.114 in the IRF 
of “P6-V11. DIFFUSE”) for front events at 30.8 GeV. 
Based on our model fit of the epoch which contains that 
highest-energy photon, we find the probability that the 
photon was associated with 3C 279 (as opposed to all 
other sources in the model including the diffuse emission 
and nearby point sources) is 88 . 6 %. 

In total, we found 10 events with estimated energies 
higher than 20 GeV within an 0°.25 radius centered at 
3C 279. All events lay within a 95% containment radius 
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Figure 3. Plots of the integrated 7-ray flux vs. photon index of 3C279 measured in week-long bins for energies above 200 MeV (A), 
between 200 MeV and 1 GeV (B) and above 1 GeV (C). Only points with TS >10 are plotted (gray points) . The blue dotted horizontal 
lines indicate average photon indices of those data for each energy band. Red points show the photon indices resulting from an analysis 
where the data were sorted in five bins using week-long fluxes for each energy band. For the red points, the horizontal bars indicate the 
ranges of the week-long flux bins while the vertical bars indicate 1 a statistical errors. The insets show enlarged views of the red points as 
well as the average photon index of each energy band. 


of the LAT PSF from 3C 279 and remain even after the 
“data clean selection” applied. The number of expected 
background events above 20 GeV within the 0°.25 radius 
at the location of 3C 279 for the 2-year observations is 
only 0.16 events. The bottom panel of Figure 1 plots the 
arrival time distribution of those 10 events. All events 
except for two were detected between MJD 54780 and 
54900 during the high activity states. No photon above 
20 GeV associated with 3C 279 has been detected after 
MJD 54914 during the 2- year observations. 

2.4. Gamma-ray spectra 

We extracted the 7 -ray spectra using data for the en- 
tire 2 -year period and following 8 sub-periods (see also 
Table 1): [A] the initial quiescent state in the 7 - ray 
band (MJD 54682 - 54728), [Bl the first "-ray flaring 
state (MJD 54789 - 54809), C] an intermediate state 
(MJD 54827 - 54877), [D] the first 5 days of the sec- 
ond 7 -ray flaring event (MJD 54880 -54885), [E] the last 
3 days of the second 7 -ray flaring event (MJD 54897 - 
54900), [F] during the isolated (first) X-ray flaring event 
(MJD 54950 - 54960; see Section 4), [G] during the sec- 
ond X-ray flaring event (MJD 55042 - 55045; see Sec- 
tion 4). and [H] a quiescent state (MJD 55240 - 55319). 
Those sub-periods were also selected taking into account 
observations in other energy bands. Spectral energy dis- 
tributions in the 7 -ray band for each sub-period are pre- 
sented in Figure 4. Each 7 -ray spectrum was modeled 
using a simple power law (PL; dN/dE oc E~ r ), a bro- 
ken power law (BPL; dN/dE oc E~ Pl for E < Z?brk and 
dN/dE oc E~ 1 ' 2 otherwise), and a log parabola (LogP; 
dN/dE oc (E/Eo)~ a -P los ( E / E °)) model. In the case of 
LogP model, the parameter /3 represents the curvature 
around the peak. We note that the choice of the refer- 
ence energy Eo in the LogP model does not affect the 
determination of the other two model parameters, and 
hence we fixed it at 300 MeV. 

The best-fit parameters calculated by the fitting pro- 
cedure are summarized in Table 1. The integral fluxes 
above 100 MeV 80 derived using each spectral model are 
also included. The averaged 7 -ray spectral shape for 



Figure 4. Gamma-ray spectral energy distributions of 3C 279 
of each period as defined in the text or Table 1. 2-year averaged 
[black: filled circles], Period A (magenta: open circles), Period B 
(browmfilled triangles), Period C (orange: open triangles), Period 
D (red: filled squares), Period E (green: open squares), Period F 
(pink: filled diamonds), Period G (cyan: open diamonds) and Period 
H (blue: crosses). The vertical bars indicate 1 <r statistical errors 
while the horizontal bars indicate energy ranges of each bin. The 
best-fit spectral models are plotted as dotted lines for each period 
and their parameters are summarized in Table 1. We use the broken 
power law model for the spectra of 2-year, Period B and Period G 
because significant improvements in the spectral fits can be seen 
compared to the simple power law model (see Table 1) while the 
simple power law model is used for other periods. The lower panel 
shows the residuals, plotted as x(= (data — model)/data error) 
from the besMit models. A ‘dip’ feature at ~ 1-2 Gel in the 
spectrum of Period G (the 3rd point in cyan) is a < 2 er effect from 
the best-fit model, thus, not statistically significant. 


60 Although we use photon data from 200 MeV, the integral 
fluxes are extrapolated down to 100 MeV, which is convenient to 
compare with other 7-ray results. 
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Table 1 

Results of spectral fitting in the 7 -ray band measured by Fermi- LAT. 


Period 

(MJD) 

fitting model 3 

Gamma-ray spectrum (Ferrai-LAT) 
r/o/r ; p/r 2 u brk (Gev) 

TS 

— 2AL b 

Flux (> 100 MeV) 
(10 -7 ph cm -2 s -1 ) 

2 years 

2008 Aug 4 - 2010 Aug 6 
(54682 - 55414] 

PL ' 
LogP 
BPL 

2.38 ±0.02 
2.18 ±0.03 
2.31 ±0.02 

0.08 ± 0.01 
2.95 ±0.12 

3.5 ±0.3 

20272 

20267 

20286 

46.5 

43.0 

6.10 ±0.13 
5.18 ±0.16 
5.76 ±0.15 

Period A 

PL 

2.30 ±0.07 



757 


3.7 ±0.4 

2C08 Aug 4 - 2008 Sep 19 

LogP 

2.19 ±0.15 

0.04 ± 0.05 


797 

0.7 

3.3 ±0.5 

(54682 - 54728) 

BPL 

2.21 ±0.09 

2.82 ±0.40 

3.4 ± 0.8 

798 

2.5 

3.4 ±0.4 

Period B 

PL 

2.28 ±0.04 



3209 


19.0 ± 1.1 

2C08 Nov 19 - 2008 Dec 9 

LogP 

1.95 ±0.10 

0.13 ±0.04 


3214 

13.6 

15.0 ± 1.3 

(54789 - 54809) 

BPL 

2.00 ±0.10 

2.61 ±0.11 

1.0 ±0.2 

3215 

13.7 

15.8 ±1.3 

Period C 

PL 

2.25 ±0.04 



4107 


10.0 ±0.5 

2008 Dec 27 - 2009 Feb 15 

LogP 

2.05 ±0.08 

0.08 ± 0.03 


4110 

8.2 

8.6 ±0.6 

(54827 - 54877) 

BPL 

2.07 ± 0.08 

2.43 ± 0.08 

1.0 ±0.2 . 

4109 

8.2 

8.9 ± 0.6 

Period D 

PL 

2.36 ±0.08 



1236 


23.6 ± 2.3 

2009 Feb 18 - 2009 Feb 23 

LogP 

2.16 ±0.16 

0.09 ±0.06 


1234 

2.0 

20.2 ± 2.8 

(54880 - 54885) 

BPL 

2.25 ±0.12 

2.91 ±0.61 

2.3 ±2.1 

1235 

2.9 

21.8 ±2.6 

Period E 

PL 

2.64 ±0.32 



61 


6.3 ±2.5 

2009 Mar 7 - 2009 Mar 10 
(54897 - 54900) 

LogP 

2.64 ±0.32 

0.00 ± 0.00 


61 

0.0 

6.3 ±2.5 

Period F 

PL 

2.54 ± 0.24 



85 


3.5 ± f.2 

2009 Apr 29 - 2009 May 9 
(54950 - 54960) 

LogP 

2.54 ±0.24 

0.00 ±0.00 


85 

0.0 

3.5 ±1.2 

Period G 

PL 

2.44 ±0.13 



460 


18.8 ±2.9 

2009 Jul 30 - 2009 Aug 2 
(55042 - 55045) 

LogP 

2.37 ± 0.25 

0.03 ±0.10 


460 

0.1 

17.7 ±4.0 

Period H 

PL 

2.83 ±0.11 



sag - 


577 ±0.5 

2010 Feb 13 - 2010 May 3 

LogP 

2.56 ±0.23 

0.16 ±0.13 


399 

1.9 

2.9 ±0.6 

(55240 - 55319) 

BPL 

2.72 ±0.43 

3.47 ±0.57 

1.6 ±0.5 

399 

1.6 

3.4 ±1.7 


a PL: power law model, LogP: log parabola model, BPL: broken power law model. See definitions in the text. 
b LL represents the difference of the logarithm of the likelihood of the fit with respect to a single power law fit. 


the 2-year observation significantly deviates from a sin- 
gle power law. A LogP model is favored to describe 
the 7-ray spectral shape over the simple PL model with 
the difference of the logarithm of the likelihood fits 61 
— 2A L = 46.5 (corresponding to a significance level of 
~ 7 o) 52 , and a BPL lit yields — 2A L = 43.0. Even 
in some individual periods as defined above, the spec- 
tra deviate from a single power law: for example, the 
spectrum in the Period B yields — 2A L = 13.6. This 
is consistent with our finding in Section 2.2 that the 
spectrum above 1 GeV is significantly softer than the 
spectrum below 1 GeV. We thus conclude that the 7-ray 
spectrum significantly deviates from a simple power law. 
The spectral break in 3C 279 is not as pronounced as that 
seen in the spectra e.g., of 3C 454.3 (Ackermann et al. 
2010). One the other hand, the BPL model returns 
break energies within a few GeV range regardless of the 
flux levels as observed in other bright FSRQs, such as 
3C 454.3 and 4C+21.35 (Tanaka et al. 2011). Such a 
spectral feature could be due to 7 — 7 absorption to 
pair production by He II Lyman recombination contin- 
uum UV photons from the emission line region (see, e.g., 
Poutanen & Stern 2010), or a break in the electron dis- 
tribution (Abdo et al. 2009b). We consider the 7-ray 
emission region to be located significantly beyond the 
broad emission region (see the Discussion in Section 5), 
and this implies that the break in the electron energy 
distribution is the more likely explanation. 

61 — 2AL = — 2 log(L0/Ll), where LO and LI are the maxi- 
mum likelihood estimated for the null and alternative hypothesis, 
respectively. 

62 Because the LogP model has one more free parameter than 
the PL model has, the — 2A L distribution follows a x 2 distribution 
with 1 degree of freedom. 


3. MULTI-WAVELENGTH OBSERVATIONS AND DATA 
REDUCTION 

3.1. X-ray and Hard X-ray: Suzaku 

The Suzaku X-ray satellite (Mitsuda et al. 2007) ob- 
served 3C 279 as a part of multi-band studies of the 
object. The observations took place in two segments, 
with an interruption lasting roughly 1.5 days: (1) be- 
tween 2009 January 19, 23:19:00 and 2009 January 22, 
22:32:00 UTC (sequence number 703049010), and (2) be- 
tween 2009 January 23, 20:45:00 and 2009 January 25, 
03:00:00 UTC (sequence number 703049020). “Period 
C” (see Table 1) includes both Suzaku observations. The 
goals of the Suzaku observations were to monitor the soft- 
medium X-ray flux (0.3-12 keV) of the source with the 
X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) 
and to take advantage of the data from the Hard X- 
ray Detector (HXD; Takahaslii et al. 2007). The HXD 
consists of PIN silicon diodes for the lower energy band 
(10-70 keV) and GSO scintillators for the higher energy 
band (40-600 keV), to extend the spectral bandpass be- 
yond the energies accessible with imaging instruments 
(> lOkeV). The HXD nominal position was used for the 
observations to maximize its effective area. In the follow- 
ing analysis, the HXD/ GSO data were not used because 
there was no significant detection of the source. 

Although the observation conditions were nominal, the 
XIS1 data suffered from somewhat high and variable 
background, resulting in the total apparent counting rate 
ranging from 1 to 3 counts s -1 in source-free regions for 
the entire chip. Still, the background-subtracted spec- 
trum determined from the XISl data below 8 keV was 
entirely consistent with that from XIS0 and XIS3 and 
thus we included the background-subtracted XISl data 
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in the spectral fitting. The total duration of good data 
accumulated by the XIS instruments was 191 ks. We 
used the standard ftools data reduction package, pro- 
vided by the Suzaku Science Operations Center, with 
the calibration files included in CALDB ver. 4.3.1. For 
the analysis of spectra and light curves, we extracted 
the counts from a region corresponding to a circle with 
260"radius, centered on the X-ray centroid; we used a 
region of a comparable size from the same chip to ex- 
tract the background counts. The net count rates were 
0.47, 0.63, and 0.56 count s -1 for XISO, XIS1, and XIS3, 
respectively, with the typical count rate uncertainty in 
the entire observations of ~ 3 %. The data indicate no 
significant variability during the Suzaku observations. 

The source was also detected in the HXD/PIN data, al- 
though the signal was relatively weak. We used the stan- 
dard cleaned events, processed using the standard crite- 
ria applicable to the rev. 2.13 of the Suzaku HXD data 
processing software. This yielded 95.4 ks of good data, 
with a total count rate of 0.3 count s -1 . For the back- 
ground subtraction, we used the standard background 
files provided by the Suzaku team through HEASARC. 
We applied the standard tool hxdpinxbpi which ac- 
counts for the particle background as well as for the con- 
tribution of the Cosmic X-ray Background as appropri- 
ate for the effective area and the solid angle of the HXD. 
The net counting rate was 0.02 count s , with the for- 
mal statistical uncertainty of ~ 10%. We note that this 
formal uncertainty is probably lower than the standard 
systematic error due to the background subtraction of 
3% of the average background (corresponding to 0.01 
count s -1 ). Nonetheless, even if the additional uncer- 
tainty is included, the source was still detected by the 
HXD/PIN. 

For the spectral analysis, we used the XSPEC spectral 
analysis software. For the spectral fitting of the XIS 
data, we used the standard redistribution files and mir- 
ror effective areas generated with Suzafai-specific tools 
xisrmfgen and xissimarf gen. In the spectral fits, we 
used the counts corresponding to the energy range of 
0.5-10.0 keV for XISO and XIS3, and 0.5-8.0 keV for the 
XIS1. We used all three XIS detectors simultaneously, 
but allowed for a small (a few %) variation of normaliza- 
tion. For the HXD/PIN data, we considered the data 
in the range of 20-50 keV and used the response file 
aeJixd-pinhxnome5_20080716 . rsp. 

The source spectrum was modeled as an absorbed 
power law, with the cross-sections and elemental abun- 
dances as given in Morrison & McCammon (1983); other 
absorption models give similar results. The best-fit ab- 
sorbing column was (3.1 ± 0.5) x 10 20 cm -2 , and the 
photon index was 1.76±0.01. Inclusion of the HXD/PIN 
data in the fit did not change the fit parameters percep- 
tibly. The x 2 for the fit including the three XIS de- 
tectors and the HXD/PIN was acceptable, with 5061 
for 5023 PHA bins. The absorption inferred from the 
simple absorbed power law model is marginally greater 
than the value inferred from the radio measurements of 
the column density of the material in the Galaxy of 
2.0 x 10 20 cm -2 (with an estimated error of ~ 10%) 
(Kalberla et al. 2005). We deem the difference not sig- 
nificant, since at such small column densities, it can be 
accounted for by even small systematic uncertainty in 


the knowledge of the effective area of the XIS instru- 
ments at the lowest end of the XIS bandpass. Further- 
more, a modest additional column density is expected 
in the host galaxy of 3C279. The observed model 2-10 
keV flux is 8.0 x 10 12 erg cm -2 s _1 , with the statistical 
error of < 2%, which is probably smaller than the sys- 
tematic error resulting form the calibration uncertainty 
of the Suzaku instruments. We plot the Suzaku 3C 279 
spectra in Figure 5. 



Figure 5. Count spectra of 3C 279 measured by Suzaku XISO 
(black), XlSl(red), XIS3 (green) and HXD/PIN (blue). The model 
plotted with the data is a broken power law obtained by a fitting 
these three XISs and HXD/PIN data. The lower panel shows the 
residuals for this broken power law model. 

3.2. X-ray: XMM-Newton 

XMM-IVewton observed 3C 279 once starting on 2009 
January 21, 17:28 UT. The observation was largely de- 
void of flares (except for the period close to the end of the 
observation), and the total length of good data accumu- 
lated in the pointing was 16.8 ks. We used the standard 
Scientific Analysis System (SAS) data reduction pack- 
age, provided by the XMM- Newton Science Operations 
Center. Since 3C 279 is a relatively bright source, we 
considered only the pn-CCD data. We note here that 
the spectra and light curves taken by MOS-CCDs were 
entirely consistent with the results inferred from the pn- 
CCD data. 

For the analysis of spectra and light curves, we ex- 
tracted the counts from within 40" radius of the source; 
we used a region of the same size, from the same pn-CCD 
chip, to extract the background counts. The data indi- 
cate no significant variability during the XMM-IVetoton 
observation. The spectral analysis was performed us- 
ing the XSPEC v.12 spectral analysis software with the 
standard redistribution files and mirror effective areas 
included in the SAS package. We used the counts cor- 
responding to the energy range of 0.5-10.0 keV in our 
spectral fits. 

The source spectrum was first modeled as an ab- 
sorbed power law; the best-fit absorbing column was 
(2.2 ± 0.6) x 10 2 ° cm -2 , and the photon index was 
1.77 ± 0.03, with x 2 of 588 for 517 d.o.f. The result is 
consistent with the spectral analysis results of the Suzaku 
observations as described in the previous section, which 
were performed during the same period as the XMM- 
Newton observation. We also considered a broken power 
law model and found that the overall intrinsic source 
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spectrum hardens with increasing energy. The absorb- 
ing column was (3.4 ± 0.7) x lfr° cm -2 , and the low 
and high energy indices were respectively 1.83 ±0.05 and 
1.55 ± 0.2 with the break energy of 4.1 ± 0.8 keV. The 
resulting x 2 was 563, for 515 d.o.f. The broken power 
law model is statistically only marginally superior to the 
simple power law model, especially given that the absorp- 
tion inferred form the simple power law model is closer to 
the value inferred from Kalberla et al. (2005) . For either 
model, the 2-10 keV flux is 7.7 x 10~ n erg cm -2 s _1 , 
with a statistical error of 5%, which is probably smaller 
than the systematic error resulting form the calibration 
uncertainty of the XMM- IVewtcm pn-CCD. 

3.3. X-ray: RXTE-PCA 

RXTE carried out 321 observations between 2008 July 
3rd (MJD 54650) and 2010 August 12th (MJD 55420). 
Those include 52 observations based on the Cycle 12 
Guest observer (GO) program and 269 observations 
based on the Core program in Cycle 12-14. The fluxes 
resulting from the Cyclel2 GO observations have been 
already reported in Paper I. Most of the individual obser- 
vations have exposure times in a range from 1.0 to 2.5 ks. 
We analyzed the data from the Proportional Counter 
Array (PCA) following standard procedures using the 
rex script in HEASOFT v.6.9. Only signals from the 
top layer (X1L and X1R) of PCU2 were extracted for 
data analysis. The data were screened with the follow- 
ing data selection: source elevation above the horizon 
> 10°, pointing offset smaller than 0.02°, at least 30 
minutes away from a SAA passage and electron contam- 
ination smaller than 0.1. The background was estimated 
with standard procedures, and the detector response ma- 
trices were extracted with the RXTE tools (command 
PCARSO v.1'1.7.1). For the spectral analysis we re- 
binned the spectra into 11 channels. The spectra from 
the channels corresponding to nominal energies of 2.6 
to 10.5 keV are adequately fitted by a single power law 
model, absorbed by a fixed Galactic column density of 

2.2 x 10 20 cm -2 using the XSPEC v.12 software package. 
The value of the column density is based on our XMM- 
Newton results (in Section 3.2) and is also consistent with 
the value based on Kalberla et al. (2005). 

3.4. X-ray: Swift-XRT 

In the HE AS ARC data base 63 , there are 80 pub- 
licly available Swift X-Ray Telescope (XRT) observa- 
tions between 2008 July 3 (MJD 54650) and 2010 Au- 
gust 12 (MJD 55420), which include 32 pointings based 
on an approved GI proposal in Cycle-4 (Proposal num- 
ber: 5080069). The results of the flux history based 
on the data until 2009 May 31 have already reported 
in Paper I. Effective exposure times of these observa- 
tions range between 1 and 3 ks, but some have longer 
exposure times, for example, 8.9 ks for ID:35019007 
(MJD 54795), 22.5 ks for ID:35019009 (MJD 54797), 

20.2 ks for ID:35019010 (MJD 54799) and 15.4 ks for 
ID;35019011 (MJD 54800). The XRT was used in the 
photon counting mode, and no evidence of pile-up was 
found. The XRT data were reduced with the standard 
software xrtpipeline v. 0.12. 6, applying the default fil- 
tering and screening criteria (HEADAS package, v.6.10). 

63 http : //heasarc . gsfc .naaa.gov/cgi-biti/W3Bro»3e/swift .pi 


The source events were extracted from a circular region, 
20 pixels in radius, centered on the source position. Ex- 
posure maps were used to account for Point-Spread Func- 
tion losses and the presence of dead pixels/columns. The 
background was determined using data extracted from 
a circular region, 40 pixels in radius, centered on (RA, 
Dec: J2000) =(12h56m26s,-05°49'30"), where no X-ray 
sources are found. Note that the background contamina- 
tion is less than 1 % qf source flux even in the faint X-ray 
states of the source. The data were rebinned to have at 
least 25 counts per bin, and the spectral fitting was per- 
formed using the energy range between 0.3 keV and 10 
keV using XSPEC v.12. The Galactic column density is 
fixed at 2.2 x 10 2 ° cm -2 during the fittings as is the case 
in the RXTE data analysis. 

Figure 6 shows a scatter plot between photon index and 
flux in the X-ray band as measured by SwiftrXBT and 
RXTE-PCA. Generally, a “harder-when-brighter” trend 
can be seen. Only the highest-flux point measured by 
SwiftrXKT shows the photon index significantly harder 
(smaller) than 1.5. 



0 5 10 15 20 

X-ray flux (2 - 10 keV) 

Figure 6. Scatter plot of flux vs. photon index of 3C 279 in the X- 
ray band with the data taken by Swift-XKF (magenta) and RXTE- 
PCA (cyan). The horizontal dotted line represents the photon 
index value of 1.5. 


3.5. Ultra-Violet: Swift- UVOT 

The Swift Ultra-Violet/Optical Telescope (UVOT; 
Roming et al. 2005) data used in this analysis included 
all of the observations performed during the time in- 
terval MJD 54650-55420. The UVOT telescope cycled 
through each of the six optical and ultraviolet filters 
(V, B,U,W1, M2,W2). The UVOT photometric sys- 
tem is described in Poole et al. (2008). Photometry was 
computed from a 5" source region around 3C 279 us- 
ing the publicly available UVOT FTQOLS data reduction 
suite. The background region was taken from an annu- 
lus with inner and outer radii of 27". 5 and 35", respec- 
tively. Galactic absorption in the direction of 3C 279 
was adapted as given in Larionov et al. (2008), namely, 
A v = 0.093, A b = 0.123, A v = 0.147, A w i = 0.195, 
A m 2 — 0.285 and Aw? = 0.271. The measured mag- 
nitudes in each band during the 2-year observations are 
mv = 15.6 — 18.7 (75 data points), m B = 16.0 — 18.8 
(80 data points), mu = 15.1 — 18.0 (88 data points), 
mw\ = 15.3— 18.0 (84 data points), mM 2 — 15.4— 18.2 
(76 data points) and mw 2 = 15.5— 18.0 (81 data points). 
All observed data points Eire shown in Figure 7. 
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Figure 7. Light curves of all observed UV-optical-NearIR bands 
of 3C 279 in measured magnitude scale from 2008 August to 
2010 August, including W2 (Stct/fc-UVOT), A/2 (SwifthUVOT), Wl 
(SwiftrVVOT), U (Swift- WOT), B (SwifUlXOT), V (Katana, 
Swift- UVOT), R (Abastumani, Calar Alto, ST-7, GRT, MDM, 
L’Ampolla, Perkins, SLT, KVA, LT, San Pedro, St. Petersburg, 
Tijarafe), J (AZT-24, Kanata), H (AZT-24, LT) and K (AZT-24, 
Kanata) bands. 

3.6. Optical, Near-Infrared and Radio observations by 
GASP-WEBT 

The GLAST-AG1LE Support Program (GASP; 
Villataet al. 2008, 2009) is a project initially originat- 
ing from the Whole Earth Blazar Telescope 64 (WEBT) 
in 2007. It is aimed to provide long-term monitoring in 
the optical ( R band), near-IR, and mm-cm radio bands 
of 28 7-ray-loud blazars during the lifetime of the AGILE 
and Fermi 7-ray satellites. 

The observations of 3C 279 in the period considered 
in this paper were performed by the observatories listed 
in Table 2. The calibrated R-band magnitudes of the 
source were obtained through differential photometry 
with respect to the reference stars 1, 2, 3, and 5 by 
Raiteri et al. (1998). Near-IR data in the J, H, and K 
filters were acquired at Campo Imperatore and Roque 
de los Muchachos (Liverpool). When converting magni- 
tudes into flux densities, optical and near-IR data were 
corrected for Galactic reddening using Ab = 0.123 mag 
(Schlegel et al. 1998). We adapted the extinction laws 
by Cardelli et al. (1989), and the zero-mag fluxes by 
Bessell et al. (1998). 

For the observations between 2008 August and 2010 
August, the measured R-band magnitude ranged from 
14.87 to 17.81 (673 data points). The R-band data have 
the best time coverage among the IR-optical-UV bands in 
our data thanks to the participation of a number of tele- 
scopes. The emission shows strong variability and the ex- 

64 http://ww.oato.inaf . It/blazarB/webt 



Figure 8. Radio light curves of 3C 279 from 2008 August to 
2010 August measured at all observed radio frequencies: 345 GHz 
(SMA), 230 GHz (GARMA, SMA), 93GHz (CARMA), 43GHz 
(Noto), 37 GHz (Metsahovi), 22 GHz (Medicina), 15 GHz (OVRO), 

14.5 GHz (UMRAO), 8 GHz (Medicina, UMRAO) and 5 GHz 
(Medicina, UMRAO). 

cess variance (F v a r ; Eq. 1) of the source R-band flux (i.e., 
in linear scale) is 0.853 ±0.001. The near-IR magnitudes 
in the J, H and K bands were measured in ranges of 
mj = 14.91 — 15.59 (20 data points), run = 12.04—14.90 
(68 data points) and mx = 11.19 — 13.45 (20 data 
points). Those data points are shown in Figure 7. The 
radio flux densities were measured in ranges of F$. GHz = 
8.5— 12.4 Jy (109 data points), Fsghz = 9.1— 15.5 Jy (124 
data points), Fh.sgiiz = 10.3—19.4 Jy (118 data points), 
F22GHZ = 10-22 Jy (16 data points), F 37G hz = 10-20 
Jy (168 data points), F43GHZ = 10 — 22 Jy (20 data 
points), F230 GHz = 5.1 — 10.5 Jy (62 data points) and 
±345 GHz = 6.0 — 6.8 Jy (7 data points). The light curves 
of the radio flux densities in those bands are plotted in 
Figure 8. 

3.7. Optical and Near-Infrared: the Kanata telescope 

We performed the V, J and Ks-band photometry and 
polarimetry of 3C 279 using TRISPEC installed to the 

1.5 m Kanata telescope located in the Higashi-Hiroshima 
Observatory. 

TRISPEC has a CCD and two InSb arrays, enabling 
photopolarimetric observations in an optical and two 
near-IR bands simultaneously (Watanabe et al. 2005). 
We obtained 64, 42 and 17 photometric measurements 
in the V, J and Ks bands, respectively. A unit of 
the polarimetric observing sequence consisted of suc- 
cessive exposures at 4 position angles of a half-wave 
plates: 0°, 45°, 22°. 5, 67° .5. The data were reduced ac- 
cording to the standard procedures of CCD photometry. 
We measured the magnitudes of objects with the aper- 
ture photometry technique. We performed differential 
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Table 2 

List of observatories that contributed data to the campaign 


Observatory 

Detector/Telescope (diam.) 

Band 

Fermi 

Gamma my 
LAT (survey mode) 

> 200 MeV 

Suzahu 

X ray 
XIS 

0.5 -10 keV 


HXD/PIN 

15 - 50 keV 

XMM -Newton, 

PN 

0.5 -10 keV 

RXTE 

PCA 

3-10 keV 

Swift 

XRT 

0.6 - 7 keV 

Swift 

Ultm-Violet, Optical, Infrared. 
UVOT 

W2, M2,W\,U,B,V 

Spitzer 

IRS 

5-38 


IRAC 

3.6, 4.5, 5.8, 8.0 /im 
24, 70, 160 (sn 


MIPS 

Abastumani, Georgia# 

(70 cm) 

R 

Calar Alto#'* 


R 

Campo Imperatore, Italy# 

AZT-24 (110 cm) 

J,H,K 

Crimean, Ukraine# 

ST-7 (70 cm) 

R 

Goddard, USA# 

GRT 

R 

Hiroshima, Japan 

Kanata (150 cm) 

V,J,Ks, polarization (V) 

Kitt Peak, Arizona, USA# 

MDM (130 cm) 

R 

La Silla, Chile 

GROND (220 cm) 

g,r,i,z, J,H, K 

L’Ampolla# 

R 

Lowell (Perkins)# 

Perkins 

R 

Lulin, Taiwan# 

SLT (40 cm) 

R 

Roque, Canary Islands# 

KVA (35 cm) 

R, polarization (no filter) 

Roque, Canary Islands# 

LT (200 cm) 

R,H 

San Pedro Martir# 

(84 cm) 

R 

St. Petersburg, Russia# 

(40 cm) 

R 

Tijarafe# 

(35 cm) 

R 

CARMA, USA 

Jtadio 

(array) 

92.5, 227.5 GHz 

Mauna Kea, USA# 

SMA(8 6 m) 

230, 345 GHz 

Medicina, Italy# 

(32 m) 

5, 8, 22 GHz 

Metsahovi, Finland# 

(14m) 

37 GHz 

Noto, Italy# 

(32 m) 

43 GHz 

Owens Valley, USA 

OVRO (40 m) 

15 GHz 

UMRAO, USA# 

(26 m) 

5, 8, 14.5 GHz 


# GASP-WEBT 

Calar Alto data was acquired as part of the MAPCAT 
project : http : / /www. iaa.es / iagudo/ research /M APC AT 


photometry with a comparison star taken in the same 
frame of 3C 279. Its position is R.A.=12h56ml6.90s, 
Dec=-05°50'43.0"(J2000) and its magnitudes are V = 
13.660, J= 12.377 and Ks =11.974 (Raiteri et al. 1998; 
Cutri et al. 2003). The photometric data have been cor- 
rected for the Galactic extinction with Ay = 0.093, 
Aj = 0.026, and Ak 3 = 0.010. The measured op- 
tical and near-ER magnitudes by Kanata in the V, J 
and Ks bands during the 2-year observations spanned 
my = 15.54—17.27 (56 data points), mj = 13.00 — 14.58 
(37 data points) and m^s = U-21 — 11.47 (17 data 
points). Those data points are also shown in Figure 7. 

We confirmed that the instrumental polarization was 
smaller than 0.1% in the V band using the observations 
of unpolarized standard stars. Hence, we did not 
applv any corrections for it. The zero point of the 
polarization angle is corrected as standard system (mea- 
sured from north to east) by observing the polarized 
stars, HD19820 and HD25443 (Wolff et al. 1996). The 
polarization shows clear variability and the degree of 
polarization was measured in the range of 3 - 36 % dur- 
ing our 2-year observational campaign. As we reported 
in Paper I, we found a rotation of the polarization angle 


by 208° together with a sharp drop of the degree of 
polarization from ~ 30% down to a few %. The event 
was coincident with a 7-ray flare (Period D-E). In the 
second half of the 2-year observations, the source was 
generally in a quiet state in the optical band, and the 
degree of polarization was also relatively low. 

3.8. Optical and Near- Infrared: GROND 

The Gamma-Ray burst Optical/Near-Infrared De- 
tector (GROND; Greiner et al. 2008) mounted at the 
MPI/ESO 2.2 m telescope at LaSilla observatory in Chile 
observed the field of 3C 279 in two nights of 2008 July 
(2008 July 30 and 2008 July 31) and four nights in 2009 
January (2009 January 19 to 2009 January 22). In each 
observation, a total of 4 images in each g'r'i'z' filter with 
integrations times of 35 s and 24 images of 10 s exposure 
in each JHK S were obtained simultaneously. 

GROND optical and near-IR data were reduced in 
standard manner using pyraf/IRAF (Tody 1993) simi- 
lar to the procedure outlined in Kruhler et al. (2008). 
The stacked images of each observation were flux cal- 
ibrated against GROND observations . of SDSS fields 
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(Abazajian et al. 2009) taken immediately before or af- 
ter the field of 3C 279 for the optical g'r'i'z', and mag- 
nitudes of 2MASS field stars (Skrutskie et al. 2006) for 
the JHK S filters. All data were corrected for the ex- 
pected Galactic foreground reddening of E^b-v) — 0 029 
according to Schlegel et al. (1998). Results of GROND 
observations are summarized in Table 3. 

3.9. Infrared: Spitzer 

We observed 3C 279 with Spitzer Infrared Spectro- 
graph (IRS), Multiband Imaging Photometer for Spitzer 
(MIPS) and Infrared Array Camera (IRAC) at several 
epochs in 2008 and 2009 under the Spitzer program 
PID50231 (PI A. Wehrle; see Table 4). The observations 
were conducted once with each instrument in 2008 July 
and August and approximately daily during the instru- 
ment campaigns in the 2009 February-March visibility 
window. 

For the IRS observations, high-accuracy blue peakup 
observations on a nearby star were used to center the 
spectrograph slit on the target. 3C 279 was observed 
with the low- resolution SL2, SL1, LL2, and LL1 modules, 
for 3 cycles of 14 seconds at each of two nod positions. 
Spitzer - IRS data reductions began with S18.7 Spitzer Sci- 
ence Center pipeline-processed, background subtracted 
data. The background was removed by subtracting the 
alternate nod for each pointing. Additional processing 
steps were applied to clean bad data, remove fringes, and 
match and trim spectral orders. First, we cleaned bad, 
rogue pixels using the Spitzer Science Center procedure 
IRSCLEAN V2.0. One- dimensional spectra were then 
extracted using the standard point-source aperture and 
flux calibration in SPICE ver. 2.3. We used a custom 
spectral defringing tool to remove fringes introduced by 
the pointing-dependent instrumental flatfield. This tool 
uses a predetermined flatfield fringing correction func- 
tion, which is shifted to match and remove the observed 
fringes in the spectrum. Spectral orders were trimmed, 
and the SL2 and SL1 orders were scaled up by a factor of 
1.06 to empirically correct for pointing-dependent point- 
source slit losses. Finally, the nod-spectra were averaged 
and combined into a single spectrum covering 5.2-35 pm 
rest wavelength. Figure 9 shows the reduced IRS spectra 
in the sF v representation. 

We used the pipeline MIPS images (ver. 18) for aper- 
ture photometry using 13", 35" and 50" radius for 24, 70 
and 160 pm bands, respectively, with aperture correc- 
tions from Tables 3.13, 3.14 and 3.16 of the MIPS Data 
Handbook (ver. 3.2). No 160 pm data were obtained 
in August 2008 because the array was not cooled dur- 
ing that campaign. 3C 279 has very low ecliptic latitude 
(0.2°), hence, the observed transients can be attributed 
to passing asteroids which appeared in various MIPS im- 
ages. We used 6" radius apertures for IRAC photome- 
try on the pipeline data (ver. 18) with aperture correc- 
tions tabulated in Table 5.7 of the IRAC Data Handbook 
(ver. 3.0). 

The 5pifzer-MIPS photometric repeatability and abso- 
lute calibration uncertainties at 24 pm are respectively, 
0.4% and 4%; at 70 pm, 4.5% and 5%; and at 160 pm, 
5% and 12% (Engelbracht et al. 2007; Gordon et al. 
2007; Stansberry et al. 2007). We therefore adopt over- 
all uncertainties of 10%, 10% and 20% at 24, 70 and 
160 /im, respectively. No color correction has been ap- 



Figure 9. Spectral energy distribution of 3C 279 in the infrared 
band measured by Spitzer- IRS. IRS spectra from highest to lowest 
are on 2009 March 3, 2009 March 4, 2009 March 5, 2009 March 6, 

2009 March 7, 2009 March 9, and 2008 August 16. Representative 
flux data points measured by Spitzer-MIPS (circles) and Spitzer- 
IRAC (triangles) are also included. 

plied to the data because the slope and smoothness of the 
spectrum over the bandpasses are not known. Figure 10 
describes the MIPS flux history during the 6 epochs from 
15 to 20 February 2009 together with R-band flux for 
comparison. No significant flux variation is found in all 
MIPS bands during those epochs, which include period 
D. 

The Spitzer-IRAG calibration uncertainty is 3 % overall 
and has photometric repeatability of 1.5% (Reach et al. 

2005) . We adopt the overall IRAC calibration uncer- 
tainty of 3%, but note the following characteristics of 
our images. In our IRAC frames, two standard compar- 
ison stars used in blazar monitoring were visible in the 
3.6pm images (Star 1 and Star 2) 65 . One comparison 
star, Star 2, was visible in the 4.5, 5.8, and 8 pm images, 
located at the interstice of the chopping regions where 
the data are noisier than elsewhere. The spacecraft ori- 
entation, and hence the chopping orientation, was 180 
degrees different between 2008 July-August and 2009 
March. The standard deviations in comparison Star 2’s 
measurements in 2009 March at 3.6, 4.5, 5.8 and 8 pm 
are 0.68 mJy, 0.08 mJ’-, 0.08mJy and 0.09 mjy (5%, 1 %, 

2 % and 3 %), respectively. The high 3.6 pm standard de- 
viation was affected by a single high value on 11 March 
2009, for which we found no obvious cause; excluding 
that value resulted in a standard deviation of 0.06 mJy 
(0.5%). In contrast, the flux of 3C 279 shows a steady 
decrease of 10%, 12%, 14% and 13% at 3.6, 4.5, 5.8 
and 8.0 pm, respectively, during the 6 epochs from 2009 
March 10-16 as shown in Figure 11. 

3.10. Radio: CARMA 

Observations were obtained at mean frequencies of 92.5 
and 227.5 GHz using the Combined Array for Research 
in Millimeter-wave Astronomy (CARMA; Bock et al. 

2006) . In all cases the nominal signal to noise ratio 
exceeded 400 and calibration uncertainties dominated 

65 see http : //vvw. lsw . uni- heidelberg . de /pro j ect s/ extragalact ic/charts, 
Raiteri et al. (1998), Villata et al. (1997) and 
http: //quasar . colgate.edu/- tbalonek/optical/3C279coopstars . gif 
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Tfeble 3 

Results of GROND observations. 



2008 July 31 

2009 .January 19-22 


filter 

AB magnitude 

flux |mJyJ 

AB magnitude 

flux |mJy] 

A/A(V)* 

sf 

17.62 ± 0.05 

0.324 ±0.015 

16.06 ±0.05 

1.37 ±0.06 

1.23 

r' 

17.24 ±0.05 

0.462 ±0.021 

15.49 ±0.05 

2.31 ±0.11 

0.80 

i' 

16.83 ±0.05 

0.671 ±0.031 

15.05 ± 0.05 

3.47 ±0.16 

0.62 

z' 

16.67 ± 0.05 

0.776 ± 0.036 

14.80 ±0.05 

4.37 ± 0.20 

0.45 

J 

16.05 ± 0.06 

1.387 ±0.077 

14.16 ±0.06 

7.86 ± 0.44 

0.29 

H 

15.60 ±0.07 

2.098 ±0.133 

13.62 ± 0.07 

12.90 ±0.81 

0.18 

K 

15.18 ±0.09 

3.062 ± 0.249 

13.19 ±0.09 

19.18 ± 1.56 

0.14 


Note. — Results of both AB magnitude and flux are corrected for Galactic extinction. 
No significant daily variability was observed during the observations in 2009 January. 
a dereddening factors for correction of Galactic extinction. 


Table 4 

Spitzer observation log. 


Instrument 

Start time (UTC) 

(MJD) 

Duration (min) 

ObsID 

MIPS 

2008.7.31 11:09:06.3 

54678.4647 

9.19 ■ 

27434240 


2009.2.15 07:08:08.6 

54877.2973 

14.61 

27438592 


2009.2.16 23:12:32.0 

54878.9670 

14.61 

27438080 


2009.2.17 20:42:36.5 

54879.8629 

. 14.61 

27438848 


2009.2.18 19:39:30.8 

54880.8191 

14.60 

27438336 


2009.2.19 14:19:04.4 

54881.5966 

14.60 

27439360 


2009.2.20 09:23:37.4 

54882.3914 

14.60 

27439104 

IRS 

2008.8.16 14:32:42.6 

54694.6060 

TL92 

27425024 


2009.3.3 12:22:35.7 

54893.5157 

18.15 

27435776 


2009.3.4 19:53:03.5 

54894.8285 

18.15 

27437312 


2009.3.6 10:50:02.7 

54896.4514 

18.16 

27436544 


2009.3.7 00:46:03.4 

54897.0320 

18.16 

27435520 


2009.3.8 12:08:53.5 

54898.5062 

18.18 

27437056 


2009.3.9 18:39:51.2 

54899.7777 

18.18 

27436288 

mxc 

2008.8.17 06:57:26.3 

54695.2899 

RT53 

27429632 


2009.3.10 19:26:05.5 

54900.8097 

10.69 

27433216 


2009.3.11 22:05:54.7 

54901.9208 

10.70 

27432448 


2009.3.13 05:13:56.9 

54903.2180 

10.71 

27433728 


2009.3.14 04:12:41.8 

54904.1755 

10.72 

27432960 


2009.3.15 03:04:16.9 

54905.1280 

10.73 

27433984 


2009.3.16 14:58:46.6 

54906.6241 

10.75 

27433472 
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Figure 10. Light curve of 3C 279 at 24 pan, 70 pan and 160 pm 
measured by Spitzer-MIPS. The error bars correspond to 10%, 10% 
and 20%, respectively, for each band as mentioned in the text in 
Section 3.9. Optical /-{-band data taken by the ground-based tele- 
scopes are also plotted in red color for comparison. 


the errors. The source is bright enough to permit self- 
calibration on timescales of less than a minute and so 
atmospheric decorrelation was not expected to affect our 
results significantly even at the long baselines. However, 
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MJD - 54000 


Figure 11. Light curve of 3C 279 at 3.6 pm, 4.5 pm, 5.8 pm and 
8.0 pm measured by Spitzer-IRAC. The error bars correspond to 
3% for all IRAC bands as mentioned in the text in Section 3.9. 
Optical /{-band data taken by the ground-based telescopes are also 
plotted in red color for comparison. 


observations in poor weather were not used due to the 
difficulty of reliably measuring pointing offsets in these 
conditions. 

Data calibration and analysis was done with the 
MIRIAD software package (Sault et al 1995). Flux den- 
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sities were determined by first using phase self calibra- 
tion with a short enough averaging interval to avoid any 
atmospheric phase de-correlation, then the flux density 
was determined from the vector average fringe ampli- 
tude at the position of 3C 279 over all baselines. For 
a strong point source such as 3C 279 this provides very 
robust and unbiased amplitude estimate independent of 
the weather or the interferometer baselines. We rely 
on regular system temperature measurements to provide 
flux calibration relative to the the fixed system sensi- 
tivity. The absolute flux calibration of CARMA obser- 
vations is usually quoted as 10-15%. However, based 
on measurements made on the blazar 3C 454.3 we esti- 
mated the relative flux calibration at each frequency to 
be within 5% at 3 mm and 10% at 1 mm. The radio 
fluxes at 92.5 and 227.5 GHz measured by CARMA cor- 
respond to F92.5 ghz = 11.7 — 19.7 Jy (14 data points) 
and .F227.5 ghz = 6.3 — 9.2 Jy (14 data points). Figure 8 
includes the flux history of those radio data. 

3.11. Radio: OVRO 40m 

The Owens Valley Radio Observatory (OVRO) 40 m 
radio data were collected as part of an ongoing long- 
term, fast-cadence 7-ray blazar monitoring campaign, 
described in detail in Richards et al. (2011). Flux den- 
sities were measured in a 3 GHz bandwidth centered on 
15.0 GHz using dual, off-axis 2f5 FWHM beams with 
12.'9 5 separation. Dicke switching against a blank sky 
reference field to remove gain fluctuations and atmo- 
spheric and ground contamination were used. Flux den- 
sities from this program are found to have a minimum 
uncertainty of 4 mJy (mostly thermal) and a typical un- 
certainty of 3 % for brighter sources. During the period 
included here, 3C 279 was observed as a pointing cali- 
brator. The flux density scale was referred to the value 
for 3C 286 (3.44 Jy at 15 GHz; Baars et al. 1977) with a 
scale uncertainty of about 5%. The radio flux at 15 GHz 
measured by OVRO was ranging from 11.1 to 18.0 Jy 
among 124 data points during the 2-year observations. 
The light curve of the OVRO radio data is also plotted 
in Figure 8. 

4. RESULTS OF THE MULTI-' WAVELENGTH 
OBSERVATIONS 

4.1. Correlations of light curves in various bands 

The multi-band light curves of 3C 279 are presented in 
Figure 12. They include (a) 7-ray flux above 200 MeV 
(Fermir LAT), (b) 2-10 keV X-ray flux measured by Swift- 
XRT and R.XTE- PCA, (c) optical- UV fluxes in R- band 
(GASP), V-band (Svdft, -U\'OT and Kanata) and W2- 
band (Swift- UVOT), (d,e) degree and angle of optical 
polarization (Kanata and KVA) and (f) radio fluxes in 
the 230, 37, 15 and 5 GHz bands (GASP, CARMA, and 
OVRO . We note that the X-rav fluxes determined by 
Suzaku and XMM-iVeurfon are entirely consistent with 
those plotted in Figure 12. The extensive data set ob- 
tained in manv bands for 3C 279 allows us to make gen- 
eral statements regarding the relative flux variability in 
various spectral bands, and the relationship of the time 
series to each other. The first such feature of the multi- 
band light curves is a general - although not exact - trend 
where the IR through optical emission seems to be cor- 
related with the 7- ray flux. We calculated the Discrete 
Correlation Function (DCF; Edelson Krolik 1988) to 


quantify the correlation of the flux variations between 
the 7-ray and other bands, and to determine whether we 
can measure any clear lag between the bands. 

For the DCF calculations, we use the 7-ray fluxes aver- 
aged over an interval of 1 day as shown in the top panel of 
Figure 1. The resulting DCF between 7-ray and optical 
22-band fluxes is shown in Figure 13. Positive values of 
‘r’ correspond to flux variations in the 7-ray band lagging 
flux variations in the other bands. In the DCF between 
'-ray and optical R- band fluxes, a. peak can be seen close 
to zero lag. We fit the DCF data points in the range be- 
tween —30 and 5 days using a Gaussian function of the 
form DCF(t) = C max x exp[(r-ro) 2 /cr 2 ], where C rnax is 
the peak value of the DCF, to is the time at which the 
DCF peaks, and a is the Gaussian width of the DCF. 
The fit yields a position of the peak at tq — —10.7 ± 0.7 
days, corresponding to a value of C max = 1.07+0.03 with 
a dispersion of a = 19.4 ±1.4 days. The result implies 
that the optical emission is possibly delayed with respect 
to the 7-ray emission by about 10 days. 

In the framework of the one-zone synchrotron + 
external-radiation Compton (ERC) models, the same 
electron population, of roughly the same energies, is re- 
sponsible for the radiation in both the optical and 7-ray 
bands. There, the observed lag can result from different 
profiles of the decreasing magnetic and radiation energy 
densities along the jet: we show that idea quantitatively 
in Appendix A. As is shown there, a very steep drop of 
the external radiation energy density is required to ex- 
plain the lag in a conical jet with magnetic field B' oc 1/r 
where r is the distance along the jet. This condition can 
be relaxed in the scenario involving the re-confinement 
of a jet (e.g., Daly & Marscher 1988; Komissarov & Falle 
1997; Nalewajko & Sikora 2009). In such a case, the 
magnetic field intensity can drop more slowly than 1/r. 
If the lag of the optical emission is confirmed, the ap- 
plication of the results in Appendix A to the ~ 10-day 
lag may imply the location of the active “blazar zone” 
at distances of a few pc in agreement with those pos- 
tulated to explain the optical polarization swing (Pe- 
riod D-E) in terms of a region containing an enhanced 
density of ultra-relativistic electrons propagating along a 
curved trajectory (Paper I). It is worth noting that simi- 
lar 7-ray/optical lags have been reported during the out- 
bursts of 3C 279 in early 1999 Hartman et al. (2001b), of 
PKS 1502+106 in 2008 (Abdo et al. 2010c), ofPKS 1510- 
089 in early 2009 (Abdo et al. 2010e; D’Ammando et al. 
2011) and of AO 0235+164 in late 2008 (Agudo et al. 
2011a; Ackermann et al. 2012). On the other hand, no 
significant lags between 7-ray and optical signals have 
been detected in 3C 454.3 in late 2008 (Bonning et al. 
2009; Jorstad et al. 2010), in 3C 66A in 2008 Octo- 
ber (Abdo et al. 2011) and in OJ 287 in 2009 Octo- 
ber (Agudo et al. 2011b). Based on investigations of 
long-term light curves of 3C 454.3 during 2008 2010, 
Raiteri et al. (2011) have shown that the optical and 7- 
ray flux variations are not always simultaneous and have 
proposed a geometrical scenario to explain the change 
in the 7/optical flux ratio during the outburst peaks in 
3C 454.3. It is expected that the on-going multi-band 
monitoring of blazars will enable us to quantify such lags 
and find out how common they are. 

Different behavior is apparent in the radio flux, where 
the energies of radio-emitting electrons are very different 
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Figure 12. Multi-band light curves of 3C 279 for 2 yeare from 2008 August to 2010 August, (a): Gamma-ray flux above 200 MeV 
averaged over 3 days, (b): X-ray flux between 2 and 10 keV measured by Sun/t-XRT and ItXTE-VCA. (c): UV-optical fluxes in R-band 
(red), V’-band (blue) and IV 2-band (magenta), (d): Polarization degree in the optical band, (e): Polarization angle in the optical band. 
The horizontal dashed lines refer to the angle of 50° and —130°. (e): Radio fluxes in 230 GHz band (magenta), 37 GHz band (orange), 15 
GHz band (blue) and 5 GHz band (green). All X-ray, UV and optica! data are corrected for tlie Galactic absorption. 
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Figure 13. Discrete Correlation Function (DCF) derived for the 
7 -ray and optical R bands. Positive values of ‘r’ correspond to flux 
variations in the 7 -ray band lagging flux variation in the optical 
band. The red curve represents a Gaussian fit to the data between 
—30 and 5 days. See the text for the fitting results. 
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Figure 14. DCF between derived for the 7 -ray and X-ray bands. 
Positive values of V correspond to flux variations in the 7 -ray band 
lagging flux variation in the X-ray band. 


from the energies of the electrons involved in produc- 
ing the observed optical and 7-ray emission. • Variabil- 
ity appears to be much less rapid, and the excess vari- 
ance (F v a,.; see definition in Eq. 1) in the radio regime 
is quite modest: for instance, 0. 145 dt 0.004 at 37 GHz, 
0.165 ± 0.001 at 15 GHz and 0.104 ± 0.001 at 5 GHz. 
Those values are significantly less than ones in the 7- 
ray or optical bands. This suggests that the synchrotron 
emission from the 7-ray emitting region is self-absorbed 
at these wavelengths. The observed radiation is produced 
at much larger distances, where the light-travel effects 
smear out the sharp, rapid variability patterns observed 
in the optical and 7-ray bands. 

Perhaps the most surprising behavior - and difficult to 
explain in the context of simple, one-component, single- 
zone models - is the relationship of the X-ray fight curve 
to those in the IR-optical or 7-ray bands. In Paper I, we 
reported that the X-ray time series exhibits a relatively 
rapid, symmetrical flare at ~ MJD 54950 (Period F) with 
a duration of ~ 20 days, which is not accompanied with 
any prominent IR./optical or 7- ray flares. As we argued 
in Paper I, the hard (rising in vF„ representation) X-ray 
spectrum is unlikely to be the “tail” of the synchrotron 
emission, but instead, it is more likely to be produced by 
the low-energy end of the electron distribution radiating 
via inverse Compton process. 

The continuing monitoring of the object in the X-ray 
band revealed another X-ray flare at ~MJD 55040 (Pe- 
riod G), ~ 90 days after the first X-ray flare. The sep- 
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aratior of the two X-ray flares is remarkably close to 
the temporal separation of the two 7-ray flares, with the 
two pairs delayed with respect to each other by ~ 155 
days. Figure 14 presents the calculated DCF between 
7-ray and X-ray fluxes, which shows a modest peak at 
~ —155 day with a correlation coefficient of 0.6-0. 7 and 
indicates no correlation between the the 7-ray and the 
X-ray bands with zero lag. While confirming the phys- 
ical connection of the two pairs would be very impor- 
tant, we cannot currently envision any situation where 
the two would be causally connected: the 155-day lag 
would imply the distance of the X-ray flare production 
~ 155P? light days ~ 50(rj/20) 2 pc and at such a dis- 
tance should be accompanied by radio flares, which are 
not seen in our data. In such a scenario, the X-ray flares 
should be significantly broadened compared to the 7-ray 
flares, however, we observe a similar temporal structure 
in both bands. Furthermore, we note that there are some 
optical and 7-ray peaks that might well be associated 
with the second X-ray flare. Hence, it is possible that 
the two prominent 7-ray/optical flares (Period B and D), 
together with the subsequent two X-ray flares (Period 
F and G), form a sequence of 4 events separated by a 
similar time intervals. Those intervals, in turn, can be 
possibly determined by instabilities in the jet launching 
region. Here, the different broad band spectra during 
these events may result from small changes of param- 
eters, such as the jet direction, Lorentz factor, and/or 
location and geometry- of the dissipation event. 

A weak (and sporadically almost absent) correlation 
between X-rays and other spectral bands can also result 
from such processes which preferably contribute to ra- 
diation in the X-ray band. They can be related to the 
following three mechanisms/scenarios: 

1. Bulk-Compton process. This involves Compton- 
scattering of ambient optical/UV light by the cold 
(non-relativistic) electrons in the jet. This mecha- 
nism is most efficient close to the accreting black 
hole where the processes responsible for the vari- 
ability of X-rays may operate independently of 
those at larger distances and producing there vari- 
able non-thermal radiation (Begelman & Sikora 
1987). A drawback of this scenario can be that 
the bulk-Compton spectrum is predicted to have a 
similar shape as the spectrum of the external ra- 
diation field (Ackermann et al. 2012), which signif- 
icantly differs from what we observe in the X-ray 
band. 

2. Inefficient electron acceleration. Acceleration of 
the relativistic electrons at proton-mediated shocks 
is likely to proceed in two steps: in the first one low- 
energy electrons may be pre-accelerated via, for ex- 
ample, some collective processes involving protons; 
in the second step, they may participate in the first- 
order Fermi acceleration process. If under some 
conditions the electron-proton coupling is ineffi- 
cient, the fraction of electrons reaching the Fermi 
phase of acceleration will be small. In this case the 
X-rays, originating from lower energy electrons, are 
produced efficiently, while the 7-rays and optical 
radiation which involve more relativistic electrons 

are not. 


3. The X-rays can be also contributed by hadronic 
processes, specifically by the pair cascades powered 
by protons losing their energy in the photo-mesonic 
process (Mannheim & Biermann 1992). For this 
process to be efficient, it requires extreme condi- 
tions (Sikora et al. 2009; Sikora 2011), however, op- 
erating in the very compact central region, at dis- 
tances less than few hundred gravitational radii, it 
may occasionally dominate in the X-ray band. 

4.2. Broad-band, Spectral Energy Distribution 

Figures 15 and 16 show broad-band SEDs of 3C 279 
in all periods as defined in Table 1. In addition, we 
also extracted a SED using data taken on 2008 July 
31 (MJD 54678), which has a good energy coverage of 
the synchrotron emission component including Spitzer 
and GROND data, although the 7-ray data by Fermi- 
LAT are not available at that time because this was be- 
fore the beginning of normal, all-sky science observations 
with Fermi-LAT. Both SEDs for Period D (2009 Febru- 
ary: corresponding to the brightest 7-ray flare coinci- 
dent with the optical polarization swing) and Period F 
(2009 April: corresponding to the first isolated X-ray 
flare) have already been partially reported in Figure 2 
of Paper I. New 5jnfzer-MIPS data points are included 
in the SED for Period D in this paper. In Period C, 
there are observations by MAGIC, which provide upper 
limits above 100 GeV (Aleksic et al. 2011). For compar- 
ison, we also include very-high-energy 7-ray fluxes de- 
tected with the MAGIC telescope in 2006 February as 
gray points (Albert et al. 2009) in Figure 15. 

This is the richest set of time-resolved spectra ever 
collected for this source. The spectral coverage of the 
synchrotron bump is unprecedented, allowing us not only 
to constrain the parameters of the emission models, but 
also to study their time evolution. As we discussed in 
Section 2.4, the shape of the 7-ray spectrum deviate from 
a simple power law, in similarity to other FSRQ blazars. 
Strong variability, over one order of magnitude, is evident 
in nearIR/optical/UV and 7-ray bands. This contrasts 
with the moderate variability in the radio/mm and X-ray 
bands. 

Particularly interesting is the behavior of this source 
in the mid-IR band, around ~ 10 13 Hz, where signif- 
icant spectral variability is observed. In the low state 
in Period A the mid-IR spectrum is relatively soft and 
can be extended with a power law shape to the opti- 
cal/UV band. In this case, the synchrotron component 
peaks in the mm/sub-mm band (~ 10 11 — 10 12 Hz). How- 
ever, in the high state in Periods D and E, the mid-IR 
spectrum is much harder and shows a significant cur- 
vature. In Period D, there is a clear spectral break at 
~ 3 x 10 12 Hz (~ 100 /um). The spectral index between 
the 70 pm (~ 4.3 x 10 12 Hz) and 160 pm (~ 1.9 x 10 12 Hz) 
points is 070-160 = 0.35 ± 0.23, taking into account sys- 
tematic errors described in Section 3.9. The synchrotron 
peak is located in the mid-IR band, at a frequency one or- 
der of magnitude higher than in the low state. This indi- 
cates that there are two independent synchrotron bumps, 
possibly produced at different locations. The mid-IR- 
peaking component, seen only in the IR/optical/UV flar- 
ing state, is characterized by a. strong and rapid variabil- 
ity. The mm/sub-mm peaking component is more persis- 
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Figure 15. Time- resolved broad-band spectral energy distributions of 3C 279 measured in Period A-H (as defined in Table 1) and on 2008 
July 31 (MJD 54678), covered by our observational campaigns in 2008-2010. X-ray, UV-optical-nearIR data are corrected for the Galactic 
absorption. 5-digit numbers in the panel indicate MJD of the periods. For comparison, the gray open circles in the very-high-energy 7 -ray 
band represent measured spectral points by MAGIC in 2006 February (Albert et al. 2009). 


tent and dominates when the source is in the low state. 
The complex shape of the SED between the mm band 
and the 70 fim point in Period D requires a coexistence 
of these two components. A similar scenario of multiple 
synchrotron components was investigated in the case of 
3C 454.3 by Ogle et al. (2011). 

In the X-ray band, despite the smaller variability am- 
plitude, we observe some spectral changes. In particular, 
in Periods F and G, which represent the two isolated X- 
ray flares, the spectrum is very similar and harder than 
on average. Figure 15 shows that these flares are not en- 
ergetically important. If we extrapolate the X-ray spec- 
tra with power laws to the 7-ray band, we under-predict 
the observed 7-ray flux in Periods B, C, D and H. Periods 
B, C, D coincide with the high-activity 7-ray state. This 
indicates that the X-ray flux cannot originate from the 
same emission component as the 7-ray flux, at least in 
the flaring state. Because the 7 rays are correlated with 
the optical flux but not with the X-ray flux, the 7 rays 
can be related to the mid-IR-peaking synchrotron bump 
while the X-rays may correspond to the mm/sub-mm 
peaking synchrotron bump. We explore this possibility 
when modeling the SEDs at Periods A and D in Section 
5.2. 

5. MODELING THE BROAD-BAND EMISSION 

We have fitted selected SEDs with one-zone leptonic 
models described in Moderski et al. (2003), including 
synchrotron emission and self-absorption, Comptoniza- 
tion of the local synchrotron radiation [SSC component] 
and external photons [ERC component], but also includ- 
ing the opacity due to internal pair-production. The ex- 


ternal radiation includes broad emission lines (BEL) and 
infrared dust emission (IR) . Their energy densities in the 
jet co-moving frame as functions of the distance r along 
the jet are approximated by the formulae: 


/ / \ CbelFj Ld 

w BELV r ) - 3irrl EL c[l + (r/r B Bh) 0BEL ] 


(2) 


UlR r 37rrf R c[l + (t-Mr)^*] 


(3) 


where £bel = 0.1 and £ir ~ 0.1 are the fractions of the 
disk luminosity Ld — 2 x 10 45 ergs -1 reprocessed into 
emission lines and into hot dust radiation, respectively, 
n3EL = 0.1 (Ld,46) 1//2 pc and hr = 2.5(Z<d,46) 1 ^ 2 pc 
[Ld ,46 = Ld/10 46 ] are the characteristic distances where 
such reprocessing takes place, and Tj is the jet Lorentz 
factor. The external radiation fields are approximated in 
the jet co-moving frame by Maxwellian spectra peaked 
at photon energies E' bel ~ 10 eV x Pj and E' m ~ 
0.3 eV x Tj. While the radiation density profile in the 
frame external to the jet should satisfy /3 bel(ir) < 2, it 
is not applied to the profile in the jet co-moving frame. 
This is because the transformation of radiation density 
depends on the angular distribution of external radiation, 
with radiation arriving at small incidence angles to the 
jet velocity vector being actually deboosted. This can re- 
sult in a steeper profile of the radiation density in the jet 
co-moving frame. Here, we adopt Abel = 3 (Sikora et al. 
2009) and pm = 4 (see Section 5.1). We assume a coni- 
cal jet geometry; the magnetic field, assumed to be domi- 
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Figure 16. Time-resolved broad-band spectral energy distributions of 3C 279 in each period (A-H), covered by the campaigns. The data 
points are the same as ones in Fig. 15, but are plotted in a separate panel for each period. 5-digit numbers in each panel indicate MJD of 
the observing period of each broad-band spectrum. 
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nated by the toroidal component, is taken to decline with 
distance r as B' oc 1/r. Electrons are injected with a 
double-broken power law distribution Q(l) oc 7 ' p with 
p = pi for 7 < Tbri, P = P 2 for Tbri < 7 < 7 br 2 and 
p = Pi for 7 > 7 br 2 - Their evolution, including injec- 
tion at a constant rate as well as radiative and adiabatic 
cooling, is followed over a distance A r = r/ 2 , where r is 
the position at which the injection ends. The emission is 
integrated over spherical thin shells within a conical re- 
gion of opening angle 0j = 1/Ij. The observer is located 
within the jet opening cone, i.e. 0 o b s < 6 y 

We begin by modeling the SED in Period D, which is 
the highest 7 -ray state reached by the source dining our 
observational campaigns. In Paper I, we showed that 
the flare event was accompanied by an optical polariza- 
tion swing and proposed two interpretations of this event. 
The first one involved a cloud containing ultra-relativistic 
particles propagating along a curved trajectory. The du- 
ration of the polarization swing constrains the location 
of the cloud to be at a few parsecs from the central su- 
permassive black hole, in the region where external ra- 
diation is dominated by the infrared dust emission. In 
Section 5.1, we present an ERC-IR model describing the 
SED in Period D and a physically related model of the 
SED in Period E. The second interpretation of the polar- 
ization swing involved the jet precession, which allowed 
arbitrary location of the emitting region, including the 
broad-line region. In Section 5.2, we present an ERC- 
BEL model of the SED in Period D. We show that in 
this scenario the far-IR break arises due to synchrotron 
self-absorption. We also show an ERC-IR model of the 
SED in Period A, which can explain the mm/far-IR and 
X-ray emission, as well as the low-state optical and 7 -ray 
flux levels. 

We assume the scenario where the X-ray emission is 
unrelated to the flaring component, since it showed lit- 
tle variability during the correlated 7 -ray /optical flares. 
Our one-zone models of the flaring states are fitted only 
to the IR/optical/UV and 7 -ray data, treating the simul- 
taneous X-ray spectrum as only an upper limit to the 
SSC component and the ERC component from the low- 
energy electrons. The large 7 -ray /X-ray luminosity ratio 
forces us to adopt a very hard electron energy distribu- 
tion at low energies (p = 1 ), which can be alternatively 
obtained by imposing a minimum electron Lorentz factor 

'Ymin 1 . 


5.1. Propagation scenario for the emitting region 

An intrinsically spherically symmetric emitting region 
is expected to produce the observed electric polarization 
vector aligned with the projected velocity of the emitting 
region. Nalewajko ( 2010 ) presented a simple model of its 
trajectory to explain the event of simultaneous smooth 
variations of the polarization degree and angle during the 
polarization swing which has been reported in Paper I. 
This model adopts a constant jet Lorentz factor Tj = 15 
and can be used to predict the viewing angles for a given 
observation time. For Period D we estimate 0 o bs,D ~ 
1.5°, while for Period E: 0 o bs,E ~ 2.4°. Between Periods 
D and E (A< ~ 15 days) the emitting region propagates 
over a distance A r ~ Tfc(At) ~ 2.8 pc. 

In Figure 17 we show Model D1 fitted to the SED in 
Period D at r = riR and Model El fitted to the SED in 


Period E at r = tir + A r. Model parameters are listed in 
Table 5. Both models use the magnetic field scaled to the 
same value at the distance of 1 pc. In order to explain the 
difference in the luminosity ratio of the ERC component 
and the synchrotron component, which decreased by fac- 
tor ~ 4 between Model D1 and Model El, we assume a 
distribution of the co-moving IR radiation energy den- 
sity dropping steeply with distance, adopting /Sir ~ 4. 
This corresponds to a strongly stratified torus structure, 
with a significant concentration of hot dust very close to 
the sublimation radius (see, e.g., Mor & Netzer 2012). 
We should note that, although the relatively soft 7 -ray 
spectrum was observed at Period E (r = 2.64 ± 0.32, 
see Table 1), the peak of the ERC-IR component in the 
Model El falls at ~ 800 MeV, in the Fermi- L AT band. 

The far-IR spectral break in Period D requires a 
sharp break in the electron distribution function at 


bV keV MeV GeV 



v|Hz] 


Figure 17. Emission models Dl (red line) and El (blue line) fitted 
to the spectral states at Periods D and E, respectively. Periods D 
and E correspond to the first 5 days and the last 3 days of the 
7 -ray flaring event accompanied by an optical polarization change, 
respectively. Those models adopt our ‘propagation scenario’, where 
external radiation is dominated by infrared dust emission. We 
assume the X-ray emission is not related to the flaring component, 
and consider the X-ray fluxes as only upper limits to the SSC and 
the ERC components during the flaring event. See the text in 
Section 5 and 5.1 for details of the models and Table 5 for model 
parameters. 
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Figure 18. Emission Models A‘2 (green line) and D2 (red line) fit- 
ted to the spectra! states at Periods A and D, respectively. Period 
A represents a quiescent state, and Period D corresponds to the 
7 -ray flaring event accompanied by an optical polarization change. 
Those models adopt our ‘jet precession scenario’, which assumes 
the 7 -ray flaring e-rent (Period D) occurs within the broad-line 
region while the low-steady emission component (Period A) is gen- 
erated outside the broad-line region. See the text in Section 5 and 
5.2 for details of the models and Table 5 for model parameters. 
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Table 5 

Parameters of emission models. 


Model 

Dl 

El 

D2 

A2 

ext. rad.® 

"TET" 

IR 

BEL 

IR 

r [pel 

1.1 

3.9 

0.045 

3.9 

R [pc] b 

0.07 

0.26 

0.0023 

0.19 


15 

15 

20 

20 

h 

3.8 

3.8 

2.9 

2.9 

flobs [°] 

1.5 

2.4 

1.7 

1.7 

Slpc [G] c 

0.14 

0.14 

0.15 

0.15 

K.-t I 10 4 er 8 cm 3 ] d 

78 

0.97 

8 x 10 4 

1.8 

Tbri 

800 

800 

170 

440 

7br2 

6500 

5000 

1000 


Pi 

1 

1 

1 

2.2 

P2 

2.5 

2.6 

2.4 

3.4 

PS 

6 

4.2 

7 



a A dominant component of the external radiation. IR: 
infrared dust emission, BEL: broad emission lines. 
b Radius of source emission region. 
c Magnetic field intensity at the distance of 1 pc. 
d Energy density of the dominant component of the external 
radiation in the jet co-moving frame at the given distance 
r. 

Tbri = 800. As the cooling break is expected at 7c ~ 
3m e c i /(2aTRu , ext ) ~ 660 (where R is a radius of source 
emission region and u' ext is the energy density of the ex- 
ternal radiation in the jet co-moving frame), Tbn is lo- 
cated just within the fast-cooling regime. The electron 
distribution in the fast-cooling regime cannot be harder 
than p— 2, hence the resulting synchrotron spectral in- 
dex a = (p— l)/2 should be larger than 0.5. In fact, the 
mean value of the observed spectral index between 70 pm 
and 160 pm is smaller than 0.5 (0:70-160 = 0-35 ± 0.23), 
which cannot be explained if the electron cooling is effi- 
cient. However, the uncertainty of the measurement does 
not allow us to reject this scenario. 

5.2. The jet precession scenario: two synchrotron 
components 

Alternatively, if the jet precession can cause the ob- 
served 7-ray flare event with the polarization swing, the 
7-ray/optical emission can be generated much closer to 
the central black hole, even within the broad line region 
(see also in Paper I). Therefore, we also attempted to 
model Period D placing the emitting region at ubel ■ For 
Tj = 15, with model parameters fitted using the syn- 
chrotron and ERC components, the X-ray flux is over- 
produced by the SSC process. To alleviate this problem, 
we increased the jet Lorentz factor to Tj =20. In Figure 
18, we show Model D2 with parameters listed in Table 5. 
The magnetic field strength scaled to the distance of 1 pc 
is almost the same as the value in Model Dl. Because of 
a smaller size of emission region and higher energy den- 
sity of the locally produced synchrotron radiation, the 
synchrotron self-absorption is able to produce a spectral 
cut-off at a higher frequency of ~ 3 x 10 12 Hz (~ 100 pm), 
consistent with the far-IR break. This interpretation has 
an advantage that it also could explain the observed hard 
spectral index between 70 pm and 160 pm, even smaller 
than 0.5, independently of details of the electron energy 
distribution. 

The low-energy synchrotron component, dominating 
the mn/sub-mm band, must be produced in a much 
larger region, placing it far outside the broad-line re- 
gion. In Figure 18, we present Model A2, fitted to the 
SED at Period A. We kept the Lorentz factor and the 


magnetic field consistent with Model D2, but we set the 
source at the distance ~ 4 pc, the same as in Model El. 
This low state model of Period A can reproduce both ob- 
served X-ray and 7-ray spectra by a single broken power 
law electron distribution. The 7-ray spectral index is 
consistent with the IR/optical/UV spectral index. The 
synchrotron self-absorption is effective at ~ 10 11 Hz and 
the spectral peak is located in the mm/sub-mm band. 

Those results suggest the existence of two synchrotron 
components: one peaking in the mm/sub-mm band and 
the other peaking in the mid-IR band. The component 
with the peak in the mid-IR band is more variable, and 
can be produced at shorter distances, within the broad- 
line region, where the far-IR break can be explained by 
the synchrotron self-absorption. 

6. CONCLUSIONS 

This paper reports details of the multi-band campaigns 
on the well-known blazar 3C 279 during the first two 
years of the Fermi mission between 2008 and 2010. Some 
key results were already presented in Abdo et al. (2010a, 
Paper I). Most important of them was the coincidence 
of a dramatic 7-ray/optical flare with a change in the 
optical polarization, which we interpreted as the result 
of a compact emitting region: either propagating along 
a curved relativistic jet or located at a constant distance 
in a precessing jet. In addition, we reported on an “iso- 
lated” X-rav flare, an event without a clear counterpart 
in other bands, and taking place a few months after the 
7-ray /optical flare. Here, we extended the observation 
epoch until 2010 August yielding the best coverage of 
time-resolved SEDs ever collected for 3C 279 from radio 
through high-energy 7-ray bands. Based on those data 
we arrived at several new conclusions about the structure 
and emission models of the relativistic jet in the quasar: 

• In the high-energy 7-ray band measured by Fermi- 
LAT, the source exhibited two prominent flares 
reaching as high as ~ 3 x 10 -6 photons cm -2 s -1 
above 100 MeV in the first year while it was in a 
relatively quiescent state in the second year. No 
significant correlation between flux and photon in- 
dex has been measured in similarity to other LAT 
blazars. The 2-year averaged 7-ray spectral shape 
above 200 MeV clearly deviates from a single power 
law. The broken-power law model returns a break 
energy within a few GeV range, which does not ap- 
pear to vary with the source flux. Such behavior is 
similar to that observed in other bright FSRQs. 

• The superb temporal coverage allowed us to mea- 
sure in detail the cross correlation of the 7-ray 
and optical fluxes. The optical signals appear to 
be delayed with respect to the 7-ray signals by 
~ 10 days. Such a lag can be explained in terms 
of the simple synchrotron and inverse-Compton 
model, in the scenario where a cloud containing 
ultra-relativistic electrons propagates down the jet 
through the regions where the ratio of the external 
radiation energy to the magnetic energy densities 
decrease with distance. We have verified this idea 
qualitatively (see Appendix A), but it still needs 
specific numerical modeling to be confirmed quan- 
titatively. 
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X-ray observations reveal a pair of pronounced 
flaxes separated by ~ 90 days. Those are not con- 
temporaneous with a pair of bright 7 -ray/optical 
flares - also separated by ~ 90 days - but instead, 
are delayed with respect to the 7 -ray/optical flares 
by about 155-days. Because of such a long delay, it 
seems implausible that these events are causally re- 
lated. Instead, the possible scenarios of the X-ray 
flares may involve changes of the source parameters 
such as the jet direction, Lorentz factor, and/or lo- 
cation of the dissipation event, or may require more 
’exotic’ solutions, for instance: bulk-Compton pro- 
cess; inefficient electron acceleration above a given 
energy; hadronic processes. At this stage we can- 
not discriminate among any of those scenarios. 

The spectral coverage of the infrared band with 
Spitzer enabled us to probe the detailed structure 
of the low-energy spectral bump, attributed to the 
synchrotron radiation. Significant spectral vari- 
ability, with soft/power-law spectra in the low state 
and hard/curved spectra in the high state, as well 
as the detection of a sharp far-IR spectral break 
in the high state, strongly suggest the existence of 
two synchrotron components: one peaking in the 
mm/sub-mm band and the other peaking in the 
mid-IR band. The component with a peak at the 
mid-IR band can be responsible for emission during 
7 -ray flaring states. 

We have applied our leptonic emission model for 
the SEDs during the 7 -ray flaring state with a po- 
larization change. Adopting the interpretation of 
the polarization swing involving the propagation of 
the emitting region - that suggested in Paper I - we 
can explain the evolution of the broad-band SEDs 
from Periods D to E during the 7 -ray flaring event 
by a shift of the position of the emitting region and 
a change of the viewing angle that are consistent 
with its trajectory. We used the same distribu- 
tion of magnetic fields and only slightly changed 
electron spectra, but required a rather steep strat- 
ification of the external radiation density in the 
form of thermal emission from the dusty torus. 
In this case, the far-IR spectral break requires a 
break in the electron distribution. The observed 
■Spifzer-MIPS spectral index 070-160 = 0.35 ± 0.23 
is marginally consistent with the synchrotron emis- 
sion in the fast-cooling regime. 

We also discussed the model in which the 7 -ray 
flare is generated within the broad emission line re- 
gion at sub-pc scale from the central black hole ac- 
cording to the jet precession scenario. This model 
explains the mid-IR break during the flaring state 
of Period D by synchrotron self-absorption. Here, 
we also discussed the low-state SED in Period A 
where the mm/sub-mm- band-peaking synchrotron 
component can be dominant. The model shows the 
related ERC component can explain the steady X- 
ray emission. 


The Fermi-LAT Collaboration acknowledges generous 
ongoing support from a number of agencies and insti- 
tutes that have supported both the development and the 
operation of the LAT as well as scientific data analysis. 
These include the National Aeronautics and Space Ad- 
ministration and the Department of Energy in the United 
States, the Commissariat a l’Energie Atomique and the 
Centre National de la Recherche Scientifique / Institut 
National de Physique Nucleaire et de Physique des Par- 
ticules in Prance, the Agenzia Spaziale Italiana and the 
Istituto Nazionale di Fisica Nucleare in Italy, the Min- 
istry of Education, Culture, Sports, Science and Technol- 
ogy (MEXT), High Energy Accelerator Research Organi- 
zation (KEK) and Japan Aerospace Exploration Agency 
(JAXA) in Japan, and the K. A. Wallenberg Founda- 
tion, the Swedish Research Council and the Swedish Na- 
tional Space Board in Sweden. Additional support for 
science analysis during the operations phase is gratefully 
acknowledged from the Istituto Nazionale di Astrofisica 
in Italy and the Centre National d’Etudes Spatiales in 
France. 

The Submillimeter Array is a joint project between 
the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astro- 
physics and is funded by the Smithsonian Institution 
and the Academia Sinica. St. Petersburg University team 
acknowledges support from Russian RFBR foundation 
via grant 09-02-00092. AZT-24 observations are made 
within an agreement between Pulkovo, Rome and Ter- 
amo observatories. This paper is partly based on obser- 
vations carried out at the German- Spanish Calar Alto 
Observatory, which is jointly operated by the MPIA and 
the LAA-CSIC. Acquisition of the MAPCAT data is sup- 
ported in part by MICIIN (Spain) grant and AYA2010- 
14844, and by CEIC (Andalucfa) grant P09-FQM-4784. 
The Metsahovi team acknowledges the support from the 
Academy of Finland to our observing projects (num- 
bers 212656, 210338, 121148, and others). The Medic- 
ina and Noto telescopes are operated by INAF - Isti- 
tuto di Radioastronomia. The research at Boston U. 
was supported by NASA Fermi GI grants NNX08AV65G, 
NNX08AV61G, NNX09AT99G, and NNX09AU10G, and 
NSF grant AST-0907893. The PRISM camera at Lowell 
Observatory was developed by K. Janes et al. at BU and 
Lowell Observatory, with funding from the NSF, BU, and 
Lowell Observatory. The Liverpool Telescope is operated 
on the island of La Palma by Liverpool John Moores 
University in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias, 
with funding from the UK Science and Technology Fa- 
cilities Council. The Abastumani Observatory team ac- 
knowledges financial support by the Georgian National 
Science Foundation through grant GNSF/ST09/521 4- 
320. The research at U. Michigan has been funded by 
a series of grants from NASA and from the NSF. Spe- 
cific grant numbers are NASA grants NNX09AU16G, 
NNX10AP16G, NNXUA013G, and NSF grant AST- 
0607523. Funding for the operation of UMRAO was pro- 
vided by the U. Michigan. 

M. H. is supported by the Research Fellowships of the 
Japan Society for the Promotion of Science for Young 
Scientists. This work was partially supported by the Pol- 
ish MNiSW grants N N203 301635 and N N203 386337, 
the Polish ASTRONET grant 621/E-78/SN-0068/2007, 
and the Polish NCN grant DEC-2011/01/B/ST9/04845. 



22 


Hayashida et al. 


APPENDIX 

POSSIBLE EXPLANATION OF A LAG BETWEEN THE 7 -RAY AND OPTICAL FLARES 


As we discussed in Section 4.1, the multi-band time series imply that during the flaring activity detected in 3C 27 9 
the optical emission appears delayed with respect to the 7-ray emission. In the context of radiation models adopted 
here (Section 5), the same electron population produces optical synchrotron photons and also inverse-Compton q -rays 
in the fast-cooling regime. A lag between the optical and 7-ray flares may therefore result from different profiles of 
decrease of the magnetic energy density u' B (r) and the external (target) radiation energy density u’ ext (r) with the 
distance r along the jet, convolved with a non-monotonic profile of the electron injection rate within the outflow. 

In the fast-cooling (FC) regime, the power injected into the relativistic particles F*e, inj is immediately radiated away 
and determines the total broad-band luminosity produced by the cooled electrons Ltot.FC- Assuming a strong inverse- 
Compton dominance, i.e., the observed inverse-Compton peak (7-ray) luminosity a L 7 being much larger than the 
observed synchrotron peak (optical) luminosity L opt , one has 


If , inj (r ) oc .fltot, F'c(t') — L 7 (r) -j- -fiopt (c) ~ 7.-, (r ) 


(Al) 


while the optical luminosity is 


7/opt (f) — 


<i( r ) 

<xt(0 


L y (r) cx 


ttB ( r ) 

“ext ( r ) 


P ;,inj( r ) 1 


(A2) 


where we assumed 5 = Tj being independent on the position r along the jet. Hence, it is clear that while a maximum 
of L 7 (r) is determined solely by the injection rate ie,inj( r )> a maximum of L opt (r) may in general be quite different, 
depending on particular radial profiles of P e ,inj(r) and of the ratio [ug/u^J (r). 

As a specific illustrative example, let us assume that the dissipation region propagating down the jet injects non- 
thermal energy into radiating particles at the rate being a broad Gaussian function of distance r with a maximum at 
ro and a width of ro/\/2, 

P e , in j oc exp [-(r - r 0 ) 2 /r-p] , (A3) 


and that magnetic field and external photon field energy densities scale with r as power laws with indices /?g and /3 cxt 


u' 3 (r) oc r and u cxt (r) oc r 


(A4) 


Then one can find that L 7 has a maximum at r = ro, as expected, whereas L 0 pt attains a maximum at 


rcr = J X a 4- \/l + 2(/?ext-AB)) , (A5) 

which is larger than ro as long as /3 ex t > /3 b and thus results in the optical flare lagging the 7-ray flare. This is due 
to the fact that with the magnetic energy density decreasing less rapidly in the jet comoving frame than the external 
radiation energy density, the drop in the injection rate P 0 , inj(r) between ro and r cr is compensated by the increase in 
the ratio ju' B / «' xt ] (r). 

Let us further consider the particular values of /3 b = 2 and / 3 ex t = 4 discussed in Section 5. With such, assuming 
again the electron injection rate being a broad Gaussian function of the distance r along the jet as in the example 
above, the observed time lag between the optical and 7-ray flares A tabs can be evaluated as roughly 

4, -^ = S!< (§) (ijb) da,s ■ (A6) 

It is encouraging that a 10-day lag is expected for F- ~ 15 and ro cs 3 pc, which are the bulk Lorentz factor and the 
location of the dissipation region comparable to that inferred from our ERC-IR modeling. 
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